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an external medium 
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ABSTRACT 

Important astrophysical sources, such as gamma-ray bursts (GRBs) or tidal 
disruption events, are impulsive - strongly varying with time. These outflows 
are likely highly magnetized near the central source, but their interaction with 
the external medium is not yet fully understood. Here I consider the combined 
impulsive magnetic acceleration of an initially highly magnetized shell of plasma 
and its deceleration by the external medium. I find four main dynamical regimes, 
that (for a given outflow) depend on the external density. (I) For small enough 
external densities the shell becomes kinetically dominated before it is significantly 
decelerated, thus reverting to the familiar unmagnetized "thin shell" case, which 
produces bright reverse shock emission that peaks well after the prompt GRB. 
(II) For larger external densities the shell remains highly magnetized and the 
reverse shock is strongly suppressed. It eventually transfers most of its energy 
through pdV work to the shocked external medium, whose afterglow emission 
peaks on a timescale similar to the prompt GRB duration. (Ill) For even larger 
external densities there is no initial impulsive acceleration phase. (IV) For the 
highest external densities the flow remains Newtonian. 

Subject headings: gamma-rays burst: general — magnetohydrodynamics (MHD) 
— shock waves — ISM: jets and outflows 



1. Introduction 

The composition of relativistic jets or outflows in different astrophysical sources, and 
in particular their degree of magnetization, is highly uncertain and of great interest. Pulsar 
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winds are almost certainly Poynting flux dominated near the central source, and the same 
most likely also holds for active galatic nuclei (AGN) and tidal disruption events (TDEs) of 
a star by a super-massive black hole. In AGN and TDEs, since the central accreting black 
hole is super-massive, then even close to it the Thompson optical depth may not be high 
enough for thermal acceleration by ra diation pressure - the main competition to magnetic 
acceleration - to work efficiently (e.g.. I Ghisellinil |2 1 lh . In GRBs or micro-quasars, however, 
thermal acceleration could also work (r T ^ 1 is possible, or even likely), and the dominant 
acceleration mechanism is less clear. 

One of the most important open questions about outflows that start out highly magne- 
tized near the central source is how they convert most of their initial electromagnetic energy 
to other forms, namely bulk kinetic energy or the energy in the random motions of the 
particles, which also produce the radiation we observe from these sources. Observations of 
relevant sources, such as AGN, GRBs or pulsar wind nebulae suggest that the outflow mag- 
netization is rather low at large distances from the source. This is known as the a problem, 
namely how to transform from o ^> 1 near the source to o <C 1 very far from the source, 
where the magnetization parameter a is the Poynting-to-matter energy flux ratio. 

Different approaches to this problem have been considered so far. Outflows that are 
Poynting flux dominated near the source are usually treated under ideal MHD, axi-symmetry 
and steady-state (minaly for simplicity). Under these conditions, however, it is hard to 
achieve o < 1 (or a -C 1) far from the source that would enable e fficient energy dissipation 
in internal shocks ( IKomissarov et al.ll2009t lLyubarskyll2009ll2010al ). One possible solution to 
this problem is that the magnetization remains high (a 3> 1) also at large distances from the 
source and the observed emission is powered by magnetic reconnection rather than by interna l 
shocks f lLvutikov fc Blandfordl[20oi lLvutikovi[2006[ iGiannios k SpruitJbood lGianniosll2~00~8[ ) . 
Alternatively, the non-axi-symmetric kink instability could randomize the direction of the 
magnetic field, causing it to behave more like a fluid and enhancing ma gnetic reconnection, 
which both increase the acceleration and help lower the m agnetization ( IHeinz fc Begelman 
2000l ; iDrenkhahn &: Spruitl 120021; IGiannios fc Spruitl 120061 ) . Another option that may be 
relevant for AGN and GRBs ( ILyubarskyl l2010bl ). is that if the Poynting flux dominated 
outflow has alternating fields (e.g. a striped wind) then the Kruskal-Schwarzschild instability 
(i.e. the magnetic version of the Rayleigh- Taylor instability) of the current sheets could lead 
to significant magnetic reconnection, which in turn increases the initial acceleration resulting 
in a positive feedback and self- sustained acceleration that leads to a low a. 

While most previous works have assumed a steady state (i.e. no time dependence), here 
the focus is on the effects of strong tim e dependence - impulsive outflows t hat a re initially 
highly magnetized, under ideal MHD. iGranot. Komissarov fc Spitkovskyl (120111 hereafter 
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paper I) have recently found a new impulsive magnetic acceleration mech anism for relativisti c 
outflows, which is qualitatively different from its Newtonian analog ( IContopoulosI Il995l ). 
and can lead to kinetic energy dominance and low magnetizations that allow for efficient 
dissipation in internal shocks. Paper I focused mainly on the acceleration of an initially 
highly magnetized shell of plasma into vacuum, and only briefly discussed the effects of its 
interaction with the external medium. Here I analyze in detail the effects of its interaction 
with an unmagnetized external medium whose density varies as a power-law with the distance 
from the central source. 

Most astrophysical relativistic outflow sources, such as AGN, micro-quasars or pulsar 
wind nebulae (PWN), operate more or less steadily over long periods of time. Therefore, the 
deceleration of their outflow due to its interaction with external medium becomes important 
only at very large distances from the source (at the "hot spots" near the leading edge of AGN 
or micro-quasar jetg^ and at the wind termination shock in PWN). AGN or micro-quasar 
jets occasionally produce bright flares, which likely correspond to a sudden and short lived 
large increase in their jet power (or energy output rate). If the resulting ejected shell (or 
blob) of plasma is highly magnetized then it can accelerate by the impulsive mechanism 
found in paper I. Since it would be propagating in the evacuated channel cleared by the 
preceding long lived steady outflow from the same source, the deceleration by the external 
medium would become important only well after the acceleration is over. There are, however, 
also sources that are both impulsive and short-lived, such as GRBs, TDEs or potentially 
also relativistic outflows from giant flares in soft gamma-repeaters. In such sources the 
deceleration because of the interaction with the external medium can become important 
already during the acceleration stage, and this may have important implications for our 
understanding of these sources and the interpretation of their observations. 

The deceleration of an unmagnetized uniforrrj^l relativistic shell through its interaction 
with the external medium has been studied in the context of GRBs (ISari fc Piranlll995t ISari 
19971 iKobavashi fc Sarilboool : iKobavashi fc Zhanjliooi iNakar fc Piranll2004l ) . The main re- 
sults are summarized and extended to a general power-law (with the distance from the central 
source) external density profile in § [2j The deceleration of a uniform magnetized relativis- 



1 In such jets, at relatively small distances from the source the external medium can provide lateral pressure 
support that helps in the collimation of the jet and its early collimation induced quasi-steady acceleration. 

2 A non-uniform shell of e.jecta or relativistic wind with a power-law profile have also been consid- 
ered in other works (e.g. Blandford fc McKee 19761 Sari fc Meszarosl 2000t Nakamura fc Shigevama 20061 



Nousek et al. 2006 ; Granot fc KumarT 



2006; 



Levinson 



2010h . and can result in a temporally extended phase 



of energy injection into the external (afterglow) shock. For simplicity, however, this work is restricted to the 
case of a uniform shell of ejecta. 
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tic shell by an unmagnetized external medium has also been studied (IZhang fc Kobavashi 



2005 



Giannios. Mimica &: Aloyl l2008t iMimica. Giannios fc Aloyl l2009t iMizuno et al 



2009 



Lyutikovll201ll ). However, most of the treatments so far have assumed arbitrary initial con- 
ditions just before the deceleration radius where most of the energy is transfered to the 
shocked externa l medium , whic h can result in some unrealistic outcomes (notable exceptions 
are paper I and lLevinsonll2010i ). 



This work self-consistently considers the combined impulsive magnetic acceleration and 
deceleration by a unmagnetized external medium of an initially highly magnetized shell. The 
main results for the acceleration into vacuum of such a highly magnetized shell (paper I) are 
described in § [3j The test case that was studied in detail in paper I features a magnetized 
shell initially at rest whose back end leans against a conducting wall with vacuum in front 
of it, with initial width Z , magnetic field B , rest mass density p an d magnetization 
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4:71 PqC^ 



> 1 . 



X 



The shell is crossed by a strong, self-similar rarefaction wave essentially on its light crossing 



time so that at a radius Rq ~ Zq it reaches a typical magnetization (a) 



2/3 



and Lorentz 



factor (r) 



~ O", 



1/3 



It then becomes super-fast-magnetosonic and looses causal contact with 



the wall, resulting in a much slower subsequent impulsive acceleration phase in which (T) oc 
i? 1 / 3 . Eventually it becomes kinetically dominated at the coasting radius R c ~ RqCq, and at 
larger radii it starts coasting at a constant Lorentz factor ((r) ~ <tq) and spreading radially 
while its magnetization rapidly drops with radius ((cr) ~ R c /R). 

The combined acceleration and deceleration for an expansion into an unmagnetized 
external medium with a power-law density profile is addressed in detail in § |H The test case 
from paper I is generalized by replacing the vacuum with an appropriate external medium. 
Five distinct dynamical regimes are identified, and their main properties are derived and 
discussed. In regime I the external density is sufficiently low that early on it hardly affects 
the shell, which accelerates essentially as if into vacuum (as described above) until well after 
its coasting radius R c . By the time the effects of the external medium become important 
the magnetization is already low, so that regime I effectively reverts to the unmagnetized 
thin shell case (where both the reverse shock emission and afterglow emission peak on a 
timescale longer than that of the prompt GRB emission). In regime II the external density 
is sufficiently large that it starts to strongly affect the shell during its impulsive acceleration 
phase, while it is still highly magnetized. The shell then starts to decelerate or accelerate 
more slowly until it transfers most of its energy to the shocked external medium. In regime 
II the shell is highly magnetized all the way to its deceleration radius, and therefore this 
strongly suppresses the reverse shock (which is either non-existent or very weak) and its 
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associated emission. Thus, regime II can be thought of as a highly magnetized thick shell 
case, in which no bright reverse shock emission is expected, and the afterglow emission peaks 
on a timescale comparable to that of the prompt GRB. In regime III the external density is 
high enough that from the very start it inhibits the acceleration so that there is no impulsive 
acceleration phase, and the dynamics become essentially independent of the flow composition 
(i.e. of o"o, while a scales linearly with o"o but affects only the small fraction of the total 
energy that is in kinetic form, (1 + a -1 <C 1). The observational signatures of regime 

III are very similar to those of regime II. In regime IV the external density is so high that the 
flow remains Newtonian all along. This regime might be relevant for a highly magnetized 
jet trying to bore its way out of a massive star progenitor in long duration GRBs. Finally, 
regime II* occurs only for a highly stratified external medium for which it replaces regime 
II, and where also regimes I and III all show interesting and qualitatively different behavior 
compared to smaller stratifications. 

Table [1] summarizes the main notations and definitions that are used in this work in 
order to help the reader follow the text. The new results found in this work are compared 
to previous works in § and their implications are discussed in § O 

2. Deceleration of an unmagnetized impulsive relativistic flow 

Before generalizing the dynamics to the case of a highly magnetized outflow, I begin with 
a detailed description of the deceleration of an unmagnetized shell (corresponding to a <C 1 
where a is defined in the next section), that initially coasts and propagates relativistically 
into an unmagnetized external medium with a power law density profile. 

For simplicity I assume spherical symmetry, and that the original ejecta from the GRB 
form a uniform shell of initial Lorentz factor To and initial width Ao, were a subscript '0' 
is used to denote the initial value of a quantity. Bulk Lorentz factors (denoted by T), as 
well as the radius R and width A of the shell are measured in the rest frame of the central 
source (which is also the rest frame of the external medium, and thus serves as the lab 
frame), while thermodynamic quantities like the rest-mass density p, the number density n, 
the pressure p, and the internal energy density e are measured in the local rest frame of the 
fluid. A reasonable variation in T of 5T ~ r will result in a significant radial spreading 
of the shell from the spreading radius, R s ~ A Fq, so that its (lab-frame) width evolves as 
A ~ max(Ao, R/Tq) ~ Ao max(l, R/R s ). The ambient medium is assumed to have a power 
law mass density profile, p\ = AR~ k , where for simplicity I consider only k < 3, which is also 
the parameter range of most physical interest. Of particular interest are the cases k — 0, 
which corresponds to a constant density medium like the ISM, and k = 2, which is expected 
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for the stellar wind of a massive star progenitor. 

As the shell interacts with the external medium and sweeps it up, two shocks are formed: 
a forward shock that propagates into the ambient medium and a reverse shock that goes 
back into the shell and slows it down. The shocked shell material and the shocked external 
medium are separated by a contact discontinuity. There are thus four different regions: (1) 
unperturbed external medium, (2) shocked external medium, (3) shocked shell material, and 
(4) unperturbed shell material. Quantities at each region are denoted by the appropriate 
subscript i = 1,2,3,4. We have T 4 = r , r x = 1, and since regions 2 and 3 are separated 
by a contact discontinuity, T 2 = T 3 = T and P2 = P3 = V- Together with the shock jump 
conditions between regions 3 and 4 (for the reverse shock) and between regions 2 and 1 (for 
the forward shock), the resulting set of equations (together with the equations of state in the 
different regions) can be solved to obtain T, e, p2 and p% (as well as the Lorentz factors of 
the reverse and forward shock fronts) as a function of To and the density ratio / = Pa/ Pi of 
the unperturbed shell mate rial and external m edium. There are two limits for which there is 



a simple analytic solution (ISari fc Piranlll995l ): for / ^> Tq the reverse shock is Newtonian, 
ancjfl T fa r , while for / <C Tq the reverse shock is relativistic, T 2~ 1 / 2 T^ 2 / 1//4 and the 
relative Lorentz factor between the fluid in regions 4 and 3 is r 43 m 2~ l l 2 T 1 J 2 / _1 ^ 4 , where 

f — P± _ £ = lt~ k ( 2 ) 

J ~ pi 4irAc 2 T 2 R 2 - k A (3 - k)T 2 R 2 - k A ' 1 ; 

E = 10 53 £ , 53 erg is the (isotropic equivalent) kinetic energy of the ejecta shell, and 

i /(3 _ fc ) f 2.5 x 10 18 ^ 3 no 1/3 cm {k = 0) , 



{3-k)E 
4nAc 2 



(3) 

1.8 x lO 19 ^^ 1 cm (fc = 2) , 



is the Sedov radius where the (isotropic equivalent) swept up mass equals E/c 2 . Numerical 
values are provided for the physically interesting cases of k = 0, which correspond to a 
uniform medium of number density n = uq cm -3 (A = nm p where m p is the proton mass), 
and k — 2, which corresponds to the stellar wind of a massive star progenitor, with A* = 
A/(5 x 10 11 gr cm -1 ). It is clear from Eq. (J2J) that k = 2 is a critical value below which 
/ decreases with radius and above which / increases with radius, before the shell starts 
spreading (i.e. while A w A and is independent of radius). Since k = 2 is also a physically 
interesting value, it will be discussed separately below. The case 2 < k < 3 will also be 
briefly mentioned. We shall, however, first concentrate on k < 2. 



3 more accurately T = T (l - y/e) and T 43 = 1 + 2e, where e = 21^/7/ < 1 
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For k < 2, / decreases with radius. Thus the reverse shock is initially Newtonian, and 
becomes relativistic at a radius given by /(-Rat) = r 2 ,, or .R^ ~ min(i? r , Rn,o) with 



Rn,o 



(- 



E 



\A7tAc 2 T^A q 



1/(2-*) 



/3-fc 



L(3-A;)r 4 Ac 



l/(2-fc) 



4.2xl0 16 C 1/2 ^3 /2 -o 1/2 r 2 lT 3 - 1 / 2 cm 



(4) 

where ( = (1 + <z)/3, T 2 .5 = r /10 25 , T G rb = (1 + 2:)A /c = 30T 30 s is the observed duration 
of the GRB and 



(3 - k)E 
4nAc 2 T 2 



l/(3-fe) 



5.4 x lO 16 ^ 3 ^- 173 !" 273 



2 - cm (/•• 0) . 



-,2/(3-fc) 



(5) 



o 



1.8 x lO 14 ^" 1 !^ cm (fc = 2) 



is the radius where a rest mass S/FqC 2 of the external medium is swept up. In this work 
T denotes the observed time (at which photons reach the observer), while t denotes the lab 
frame time. The observed times corresponding to Rn,o an d Rr are 

Rn,o 



T, 



N,0 



1 + z) 



0- + z) 



bcT 2 



Rr 

bcT 2 



13C 3/2 K 



V2 -1/2 
53 n 



1 2.5^30 



1/2 



27C^3 /3 n - 1/3 r-f s (* = 0) 



(6) 
(7) 



0.089^53^^8 (A; = 2) , 

where b ~ 1 — 2 and b ~ 2 reflect the typical photon arrival times from regions 3 and 2, 
respectively, and b = 2 is used to obtain the numerical values. Two additional important 
radii are the spreading radius R s ~ TqA ~ R T ~ k R k ^~Q mentioned above (where the shell 
starts spreading radially), and the radius at which the reverse shock finishes crossing the 
shell, R A ~ / 1/2 r A ~ (EA/Ac 2 ) 1 '^ ~ m&x(R r , R Aj0 ) where R Afi ~ (ikRjr*)V(*-*) _ 
(-EAo/Ac 2 ) 1 / 1 - 4- ^. It is also convenient to define the parameter 

1 V(3-*) 

(3 - k)E 



T n EE 



4nAc 2 T 2 ^- k) A 3 - k 

-,1/3 -l/3 r -8/3 



Is p-2(4-*s)/(3-fc) 

a; 



l.SC^ no / r 2 l /3 T 30 1 (* = 0) 



5.9 x lO-^sA-^T^ (k = 2) 



and^l T = T (A /A) = (l s /A)T 



-2(4-fc)/(3-fc) 







. Note that Rr = Rro and R s = R s q since i?r 



does not depend on A and i? s depends on A rather than on A. Thus, we have 



(9) 



Note that T = £ , where £ is essentially the same parameter that was defined in lSari fc Piranl (|1995f ). 
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so that the initial relative ordering of the different radii is determined by the value of To, 
while the evolution of this ordering is determined by that of T. 



The condition T > 1 can be written as A < A cr or T < T cr where 

l/(3-fc) 



(3 - k)E 



5.4 x 10 11 E^n" 1/3 r^ /3 cm (k = 0) 



2(4 -k) 

r ( 3 - fe ) 

± n 



53^1* J- 2.5 Cm 



(3 - k)E 
4TrAc 2 A 3 - k 



l/2(4-fc) 



A f 



(3-fc) 
2(4-fc) 



1.8 x 10 9 E^A-^t 



395C 3/8 4 



53 n 



1/8^,-3/8 



30 



(fc = : 

(* = 0) 



88 C 1/4 ^53 /4 ^^ 1/4 T 3 - 1/4 (A: = 2) 



(10) 



(11) 



so that this case is often referred to as a "thin" or "slow" shell. Similarly, the case T < 1 
corresponds to A > A cr or r > r cr and is referred to as a "thick" or "fast" shell. Note 
that 

V a c J " lr cr 



T f 



-2(4-fc)/(3-fc) 



(12) 



2.1. Thin Shells 

For T > 1 (a thin or slow shell) and k < 2, the initial ordering of the critical radii 
is R s < Ra,o < Rt < Rn,o an d the shell starts spreading early orjfl so that at R > R s we 
have A ~ R/Tq ~ (R/R S )A and T ~ T (i?/i? s ) _1 starts decreasing, which leads to a triple 
coincidence, R& ~ Rr ~ -Rv with T ~ 1 at that radius (see Eq. [9]). In this case the reverse 
shock is mildly relativistic during the period when most of the energy is extracted from the 
shell, near the radius i?A ~ Rr ~ Rn or the corresponding time Tp when the reverse shock 
finishes crossing the shell. At larger times or radii, most of the energy has already been 
transferred to the shocked extern al medium and the flow app roaches the adiabatic (i.e. with 
a constant energy E) self-similar iBlandford fe McKed (119761 hereafter BM76) solution. 



For k = 2, / is initially (at R < R s ) independent of radius and f /Tq = /s/TqA = T 
(= To as long as the shell does not spread significantly). Therefore, for thin shells the reverse 
shock is Newtonian with a constant shock velocity at R < R s . However, for thin shells R s is 
smaller than all other critical radii, so that the shell begins to spread early on. Therefore, 
again at R > R s we have A ~ R/Tq ~ (R/R s )Aq an d T ~ T (-R/-R s ) _1 starts decreasing 



5 If there is no significant spreading of the shell (i.e. <5ro <C Tq) then the reverse shock will cross the shell 
while it is still Newtonian, and the energy extraction would pro ceed via a semi ste ady state of Newtonian 



shocks and refraction waves traveling back and forth in the shell (pari fe Piranlll9 95). In this case we do not 
expect to have significant radiation from the original shell of ejecta during its deceleration. 
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with radius, leading to Ra ~ Rr ~ Rn with T ~ 1 at that radius, so that the reverse shock 
is mildly relativistic by the time it finishes crossing the shell, at Tp. 

For 2 < k < 3 and T > 1, the initial ordering of the critical radii is Rn,o < Rs < 
Ra,o < Rr and / initially (at R < R s ) increases with radius (and time). Hence the reverse 
shock is initially relativistic until Rn,o (Tn,o) and then becomes Newtonian. At R > R s the 
shell begins to spread and from this point on A ~ (R/R s )Aq and therefore / and T begin to 
decrease with radius (as R k ~ 3 and -R -1 , respectively). This again leads to -Ra ~ Rr ~ Rn 
with T ~ 1 at that radius, where the reverse shock finishes crossing the shell. Here Rn is the 
radius where the reverse shock becomes relativistic again, i.e. it becomes mildly relativistic 
when it finishes crossing the shell, at Tp. 



2.2. Thick Shells 

For T < 1 (a thick or fast shell) and k < 2, the initial ordering of the critical radii 
is Rn,o < Rr < Ra,o < Rs- Since R s is the largest of the critical radii, spreading is 
unimportant, and therefore A « A , Rn ~ Rn,o, and R& ~ -Ra,o- The reverse shock becomes 
relativistic before it crosses most of the shell, and therefore in this case most of the kinetic 
energy is converted to internal energy (of the shocked shell and the shocked external medium) 
at -Ra,o corresponding to an observed time T E ~ (1 + £)^A,o/ c rBM(-^A,o) ~ (1 + z)A /c ~ 
T GRB , where T BM {R) ~ (E / 'A^f^R^' 2 is the Lorentz factor of the adiabatic BM76 self- 
similar solution. Here Rr is no longer relevant since the relativistic reverse shock implies 
that T To so that the energy in the swept up external medium of rest mass M is now 
T 2 Mc 2 <C TqMc 2 , and an external medium of rest mass much larger than M /T = E/TqC 2 
(by a factor of r 2 /r 2 « 2r // 1/2 > 1, where M is the original shell's rest mass) needs to be 
shocked in order for it to reach an energy comparable to E (and this occurs only at Ra,o)- 



I now generalize the results of I S aril ( 119971 ). which are for a unifo rm external density (k 



0), to a more general power law external density (with k < 3; see also lGranot Ramirez- Ruiz 



20111) • At T < T N ,o (R < R N ,o) we have T w r , while at T Nfi <T <T E (R N>0 < R < R A 



o 



we have T w 2- 1 / 2 rJ /2 / 1/4 , which can be expressed as 

1/4 / ttt \ l/2(4-fc) 



r ~ ( _ \ D(*-2)/4 w ( ^ ] T -(2-fc)/2(4-fc) (X o\ 

Vl67rAc 2 A ; Vl67r6 2 -Mc 4 - fc A ; 2 ' 



For k = 2, r remains constant at this stage while for < 2 it decreases with time (see 
below). Since T z = T/(l + z) oc i?/r 2 oc i?( 4 - fc )/ 2 we have d.R/rfT, = [2/(4 - k)]dR/dT z and 
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the observed rate of production of internal energy in the forward shock is 

dE dE dR 8nb 3 ~ k . ^ kr[2 u-k^2-k ,n 

' obB = = IR ^ = J^k) Ac T )T * • (14) 



where dE/dR » 4vr J R 2 p 1 ( J R)c 2 r 2 ( J R) = 4vrAc 2 r 2 J R 2 - ?£ . Substituting Eq. (JT3^ into Eq. (pjj 
we see that regardless of the value of k, the luminosity of the forward shock is constant, 
Wbs = (l+z)E/T E where T E = [2(A-k)/b)(l+z)A /c = [2(4-fc) /b]T GRB « 30(4-Ar) T 30 s 
is the time when the energy in the shocked external medium becomes comparable to E. The 
Lorentz factor at this time is independent of the initial Lorentz factor Tq, 



{A _ k) k-2 E 11/2(4-*) 



2 6-fc 7r7 4 c 2 A 3-fc 



r cr . (15) 



After the time T E most of the energy is in the forward shock, which quickly approaches 
the BM76 self-similar solution, in which its Lorentz factor scales as Tbm oc _R~( 3-fe )/ 2 oc 
r -(3-fc)/2(4-fc) > which imp i ies L inti0bs oc T~ l . 

For k = 2 and thick shell (To < 1) we have R s > Ra,o so that the shell hardly 
spreads radially (A « Aq) while it is crossed by the reverse shock. This implies that 
//Tq = Is/TqA = T ~ To, i.e. the reverse shock is relativistic and its strength (or r 43 ) is 
constant with radius until it finishes crossing the shell at Ra,o (corresponding to an observed 
time T E ). Therefore, for thick shells Rn,o and Tjv,o go to zero, and the Lorentz factor of the 
shocked fluid is constant in time, T(T < T E ) = T(T E ) ~ r cr (note that this value is T ). 
At T > T E (or equivalently, R > Ra,o) the flow approaches the BM76 self-similar solution. 

For 2 < k < 3 and a thick shell (T < 1), the initial ordering of the critical radii is 
Rr < Ra,o < Rs < Rn,o and / increases with radius (and time). Therefore the reverse 
shock is relativistic until it finishes crossing the shell at Ra,o (or T E ). Again, Y is given by 
Eq. (1TBT) at T < T E (or R < Ra ), where it increases with time (and radius) at this stage, 
while at T > T E (or R > Ra.o) it is given by the BM76 self-similar solution, Tbm(-R) ~ 
(E/Ac 2 y/ 2 R (k -V/ 2 . 



3. Acceleration of a highly magnetized impulsive flow into vacuum 

This was addressed in great detail in paper I, and here I summarize the main results that 
were derived in there. Paper I has studied, under ideal MHD, the test case of a cold (with a 
negligible thermal pressure) finite shell of initial (at t = 0) width Z (occupying — 1 < x < 0) 
and magnetization o"o = Bl/A.npoc 2 ^> 1, whose back end leans against a conducting wall (at 
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ce = — io)@ and with vacuum in front of it (at x > 0), where the magnetic field is perpendicular 
to the direction of motion. A correspondence was shown in this case between the dynamical 
equations in planar and spherical geometries. A strong rarefaction wave develops at the 
vacuum interface and propagates toward the wall at the initial fast magnetosonic speed of 
the unperturbed shell, c mSj o = (3 ms ,oc, reaching the wall at t = to = /o/c ms ,o ~ h/c. For a 
cold shell the dimensionless fast-magnetosonic speed is given by /3 ms ,o = y^o/JT+ ao) and 
corresponds to a Lorentz factor of r ms = (1 — /3 2 ls0 ) -1 ^ 2 — \A + and a dimensionless 
4- velocity of M mSj0 = r mSi0 /3 ms , = o"J /2 . In our case a > 1 so that r mSj0 w u msfi = (?l /2 > 1 
and /3 mSj o ~ 1. The rarefaction wave accelerates the shell to a typical (or weighted mean over 

1 /3 

the energy in the lab frame) Lorentz factor of (r)(£o) ~ cr while the typical magnetization 
drops to (er)(io) ~ cr^ 3 - This result has a simple explanation: as long as (a) ^> 1 and most 
of the energy is in electromagnetic form, energy conservation implies that (T)(a) ~ ao', such 
very fast acceleration can occur only as long as the flow pushes against the "wall" (or static 
source), and stops when the flow looses causal contact with it, i.e. when it becomes super- 
fast-magnetosonic, (r) ~ T ms ~ (o - ) 1 ^ 2 ~ <j\ ' 2 (r) -1 / 2 , which corresponds to (r) ~ a^ 3 and 
(a) ~ o" 2 / 3 - Such a shell is broadly similar to a uniform (quasi-) spherical outflow from a 
static source that lasts a finite time, to, during which it reaches a radius Ro ~ do, Lorentz 
factor (r)(t ) ~ (Tq^ 3 and magnetization (cr)(i ) ~ cr^^ \ being quickly accelerated from T ~ 1 
and a = a near the source. 

In a spherical steady-state flow the acceleration becomes inefficient once the flow loses 
causal contact with the static source (or "wall" ) and there is no significant subseque nt acceler- 



1 /S 

ation so that (r) ~ <r also asymptotically, at very large distances from the source ( IGoldreich &: Julian 



1970I ). For a non-spherical flow collimation can result in further acceleration up to (r) 



^ /o 2/3 

cx 9, (e.g., lLyubarskyll2009l ). where 9j is the asymptotic half-opening angle of the jet (at 



which point lateral causal contact across the jet is lost, so the center of the jet cannot push 
against the ambient material; for simplicity, factors of order unity are discarded here and 
until the end of this subsection). However, for an impulsive source, which corresponds to a 
shell of finite width Iq or a outflow lasting for a finite time to lo/c, efficient subsequent 
acceleration (at t > to) does occur. This happens since the shell pushes against itself and 
significantly expands in its own rest frame, under its own magnetic pressure (while its width 
in the lab frame remains constant, A = A'/T ~ l Q , since its comoving width A' increases 
linearly with its Lorentz factor T as it accelerates) . While in the comoving frame the expan- 
sion is roughly symmetric between the back and front parts of the shell, in the lab frame 
most of the energy remains in the front part of the the shell, resulting in a constant effective 



6 Such a "wall" can be the center of a planar shell surrounded by vacuum on both sides, which splits into 
two parts going in opposite directions, with reflection symmetry about its center, which remains at rest. 
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width (A ~ l , where most of the energy resides). 

The radial expansion of the shell in its own rest frame as its accelerates results in a 
dispersion 5T ~ (r) in its Lorentz factor. This causes the shell width in the lab frame to in- 
crease as A ~ Rq + R/ (r) 2 . Ideal MHD implies that the shell's electromagnetic energy scales 
as -Eem oc 1/A. Therefore, at the radius R c where the shell doubles its initial width, half of 
the initial magnetic energy is converted into kinetic form, so that (a) = E EM /E kin ~ 1 and 
(r) ~ (T & t this radius. Therefore, R c must correspond to the coasting radius where the accel- 
eration saturates and after which the shell becomes kinetically dominated and starts coasting 
at (T) ~ o"o. This, in turn, implies that R c ~ RqcTq, which provides the scaling of (r) with 
R during the acceleration phase: rflog(r)/<ilogi? = log[(r)(i? c )/(r)(i? )]/ log(-R c /-Ro) — 
\og(a^ 3 )/ log(<7o) = 1/3, so that (r) ~ (<t R/ Rq) 1 ^ during this phase, which ends at the 
coasting time, t c ~ t (jQ, distance l c ~ ct c w / o"q or radius R c ~ ct c ~ Ro<Tq. 

At t > t c the flow becomes essentially unmagnetized (i.e. with a low magnetization, 
a < 1), its internal (magnetic) pressure becomes unimportant dynamically, and each fluid 
element within the shell coasts at a constant speed (ballistic motion). As we have seen above, 
the shell starts spreading radially significantly in the lab frame at R c , and subsequently its 
width grows linearly with R, t or x, 



where ( c = t/t c w x/l c = R/R c (where R c ~ l c ). Moreover, the growth in the width of the 
shell causes a significant drop in its magnetization: a(t > t c ) ~ t c /t. One can summarize 
this result in terms of Co = t/to ~ x/l = R/Rq (where R$ ~ Iq) or ( c , 




(16) 




K Co < Oq 





(17) 





(18) 
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4. Acceleration and Deceleration of an Impulsive High-a Relativistic Outflow 

4.1. The general framework, and a spherical self-similar solution for k = 2 

For concreteness, let us specify to a spherically symmetric flow expanding into a power- 
law external density profile, p\ = Ar~ k , where r is the spherical radial coordinate. The 
outflow is taken to be cold (with no thermal pressure), and with a high initial magnetization, 
o"o ^ 1- The original outflow remains cold as long as it is not shocked by a reverse shock. The 
shocked swept-up external medium, however, is typically heated to relativistic temperatures. 
The motion is in the radial direction ((3 = f) and the magnetic field is tangential (r • B = 0). 

It has been shown in paper I that the relevant cold (no thermal pressure) MHD equations 
for spherical and planar geometries are identical when written in terms of the normalized, 
barred variables, which can apply to both a planar and a spherical geometry, 

{(x, b, p) (planar) , 

(19) 
(r, rb, r 2 p) (spherical) , 

where b = B/(y/4jcT) is the normalized comoving magnetic field. When there is thermal 
pressure then it violates this rescaling^. There is a convenient analytic solution for the 
re levant planar th e Riemann problem with a uniform unmagnetized external medium (paper 



I; iLyutikovl 120101 ) . which has a corresponding spherical solution according to the above 
rescaling, for k = 2. This solution would be valid within the original cold magnetized 
shell, i.e. at r < _Rcd(£), where Rct> is the radius of the contact discontinuity (CD) that 
forms@ A shock propagates into the cold unmagnetized external medium, with a shock radius 
Rsh{t), which heats the material passing through it to a relativistically hot temperature. 
Therefore, in the region between the shock front and the CD, Rcr>(t) < r < Rsh{t), the 
simple self-similar solution for the planar case where this region is uniform with the same 
pressure and velocity as the CD itself, is no longer valid in the spherical case. However, 
for the spherical case with k = 2 and a constant velocity of the CD (Fcd = const) there 
is a different self-similar solution shown in Figs. 4-6 of BM76, corresponding to k — 2, 
q = (energy injection by a constant power source), and m = (the Lorentz factor has no 
explicit time dependence, and instead depends only on the value of the self-similar variable, 



7 This occurs since in the momentum equation there is a term d r p or d x p, while in spherical geometry this 
rescaling requires p = r 2 p, which would instead give r~ 2 d r r 2 p = d r p + 2p/r, i.e. a spurious extra term. 

8 In this region, for k — 2,1 derive the expressions for the density p for the planar case, p p \, and those for 
the spherical case are given by p S ph(r,t) = (r / R )~ 2 p p i(x — r,t), where p sp h{r,t = 0) = p (r/R )~ 2 is the 
initial density profile of the spherical shell. 
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X = (1 + 2r 2 h )(l — £) where £ = r/ct, so that the Lorentz factor of the shock front, r s h, or 
the CD, r CD , are constant). In our case, this unmagnetized (a = 0) solution would apply in 
the region between the CD and the shock front, RcT>(t) < r < i? s h(t), while the inner part 
(r < £ u ct < R CD (t), where £ u is introduced below) of the global solution is given by the self- 
similar solution mentioned above for the cold magnetized shell, which can be simply scaled 
from planar to spherical geometry. Please note that the shock location, x — 1> corresponds 
to £ sh = ^»1- l/2r s 2 h and that the CD location is (from Table I of BM76) xcd ~ 1-77, 
corresponding to £cd = /3cd = 1 — XcD/2r 2 h and therefore 

-Sl = XCD « 1.77 . (20) 

1 CD 

However, the Lorentz factor of the material just behind the shock front is T(x — 1) ~ r s h/ V^, 
and therefore Tcb/^(x = 1) = (2/xcd) 1 ^ 2 ~ 1-06. This shows that the Lorentz factor of the 
shocked external medium increases only by about 6% from just behind the shock front to 
the CD (and its square increases by 13%, as can be seen in Fig. 5 of BM76). Therefore, a 
uniform Lorentz factor is a reasonable approximation for this region (and I shall occasionally 
use this approximation). Moreover, the normalized width of this region is 

-Rsh — Rcd ^ Rsh — Rcd ^ Xcd — 1 _ 1 — Xcd ^ 0-435 

-Rcd -Rsh 2r 2 h 2r 2 3D 2r 2 3D 

This spherical self-similar solution for k = 2 is a very useful starting point for the current 
discussion. It will be described it in terms of the corresponding planar solution, where a 
uniform region would correspond to an r~ 2 dependence of the density or magnetic pressure in 
the spherical solution. We are interested in an initially highly magnetized flow (er ^> 1), and 
for all cases of interest (except regime IV, which is described separately in § 14. 5p the shock 
that is driven into the external medium is (at least initially) highly relativistic, the shock 
front moving with T s h = (1 — /3 h 2 h) _1 ^ 2 ^> 1. The planar Riemann problem contains 5 regions 
(see Fig. [Q: (1) at £ > £ s h, where £ = x/ct and £ = £ s h = x^{t)/ct = (3 sh at the location of 
the shock front, there is cold, unmagnetized, unperturbed uniform external medium at rest 
with rest mass density pi, (2) at £cd < £ < £sh there is a uniform^! region of shocked external 
medium, moving at T 2 = (1 — /3|) -1 ^ 2 ~ T s h/ v2 with e 2 = 3p 2 = ^\p\<? , where at £cd = $i 
there is a contact discontinuity (CD), (3) at ^ u < £ < £cd there is a uniform region (moving 
at T 3 = T 2 with e3 = P3 = (_B 3 /r 3 ) 2 /87r = p 2 ) occupied by magnetized material (originating 
from region 5, or from the original magnetized outflow in an astrophysical context) that has 



9 As discussed above, in this region there is a deviation from the simple scaling between the planar and 
spherical cases, and the BM76 solution with m = q = and k — 2 holds there in the spherical case. 
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passed through a rarefaction wave (region 4) and is accumulating between the front end of 
the rarefaction wave, at0£ u = [/3 2 - l3 ms (/3 2 )]/[l- l3 2 /3 ms (l32)], and the CD, (4) at £ rf < £ < 
is a region with a rarefaction wave described by the self-similar solution in Appendix A of 
paper I, where £ r f = — /3 ms ,o = — [o"o/(l + ^o)] 1 ^ 2 is its tail, and (5) at £ < £ rf is the original 
unperturbed uniform, cold magnetized shell at rest with rest mass density po, magnetic field 
.Bo, and magnetization cr = B 2 /4:Tipoc 2 3> 1. 

Now, let us consider such an initial shell of finite initial width l , whose back end is 
leaning against a conducting "wall" (at x = —lo). At t = /o/ c ms,o (where t rs l /c for 
o"o ^> 1) the tail of the leftward moving rarefaction wave reaches the wall and a secondary 
right-going rarefaction wave forms that decelerates the material at the back of the flow. The 
heaooof the secondary rarefaction wave is located at £*(£) = x*(t)/cb and moves to the right 
with a dimensionless speed 



where > £*) and /3 m s(£ > £*) are given by the self-similar solution for the original 
expansion (describing a leftward moving rarefaction), since the part of the flow ahead of the 
secondary (or "reflected") rarefaction wave (£ > £*) does not "know" about the existence 
of the "wall". At this stage region 5 described above no longer exists, and a new region is 
formed behind the head of the secondary (right-going) rarefaction wave. This new region 
carries a very small fraction of the total energy as long as the magnetization at its head is 
large, cr(£*) = crop* 1 where p* = / po, which implies that this rarefaction is strong and 
significantly decelerates the fluid that passes through it (see paper I for details). Therefore, 
as long as this condition holds, most of the energy and momentum in the flow, as well as 
most of the original rest mass of the magnetized shell, remain in a shell of constant width 
~ 2l between and £ s h- 

The value of £ u is determined by pressure balance at the CD. Since both the normalized 
pressure, p = b 2 /2 = r 2 b 2 /2, and the fluid velocity are constant in the range £ u < £ < £ s h = 
f3 sh (corresponding to regions 2 and 3; see Fig. [TJ, and R s h(t) ~ Rcuit) = R(t) w ct so that 
the external density can be evaluated at either of these radii, pi[R s h(t)] ~ Pi[-RcdC0L we 



10 Here (3 ms ((3) is the dimensionless fast magnetosonic speed within the rarefaction wave (region 4), at the 
point where the flow velocity is v — f3c. 

11 Note that I refer to the rightmost point in the rarefaction wave as its head. In the original rarefaction 
wave this was at the vacuum interface while for the secondary rarefaction wave this is at the interface with 
the original rarefaction wave. 



1 dx* 
c dt 



i + /3(f*)/Mf.) ' 



(22) 
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have 



atr 2 CD (R) Pl (Ry= P2 (R)= P3 (R) = ( -MW) = ( ^Y^pl(R) , (23) 



-J VsCD / \H 

1 uniform approximation 

0.571 BM76 solution , 



(24) 



where p = p/p Q is the normalized density (i.e. p/po in the planar case and pr 2 j p§R\ in the 
spherical case) and p u = p(£u)/p~o is its value at £ u , while Eq. (|24|) holds for Tcd 3> 1- 

Although the self-similar solution at r < Rcr>(t) is strictly valid only for k = 2, for 
which Tcd is constant, we shall make the approximation that it still provides a reasonable 
description of the flow for k ^ 2, in which case Tcd, 0"cd, etc., gradually evolve with time. 

Denoting the initial shell to external density ratio by fo = po/pi(Ro), Eq. ( |23l) implies 
p u (R) = ^ - — — — r CD , a u — ^ — r CD , (25) 



where <r n = cx(£ u ). For the self-similar rarefaction wave solution in region 4 (see paper I), 

(rrf ) (v^+ v^TT) 4 = J + = (v^+vW+I) 4 « 16or 2 , (26) 

where a = cr(£) = a p = &op(£)/p~o is the local value of the magnetization parameter, and the 
Riemann invariant J + approaches a value of 16<Tq for cq 3> 1. We are interested primarily 
in the relativistic part of region 4, for which T 4 3> 1 is given bj 

2(T + 1 



(2<r + 1)/(1 + 2^) ~ 2<7 a < 1 , 

(27) 

(v^+v^TT) [ 1/2 - a>1 _ 

At 1 a o"o the Lorentz factor varies significantly with o = a^p as r 4 « o"o/2cr, while for 
er <C 1 it approaches a constant value of T 4 f» 2o"o. The transition between these two regimes 
occurs at a ~ 1 for which T 4 ~ cr [though r 4 (cr = 1) « 0.343cr , and r 4 (a = 1/8) = er ]- We 



12 The result for the bulk of the rarefaction wave (where r>l and a ^> 1) can be understood considering 
a finite shell of initial width Iq and energy (per unit area) Eq = Iq(Bq/8tt)(1 + 2/<tq) s» IoBq/8tt. After the 
passage of the rarefaction wave, the shell width becomes ~ 21q, and since it is relativistic there is an electric 
field in the lab frame that is almost equal to the magnetic field so that the shell energy is E w 21q(B 2 / 4w) . 
Now, E = Eo requires B sw Bq/2, and B = B^Tp since B/Tp = const, implying V sw 1/2/5. More generally, 
r 4 - + ^ X )/2, where Sp = [(1 + 0)/(l - /3)] 1 / 2 - [(vTT^ + v^)/(vT+^ + V^)] 2 - 
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are particularly interested in when this also corresponds to the transition between regions 
4 and 3, i.e. r 4 (£ u ) = T 2 ~ o"o an d ~ 1> which according to Eq. fl25|) corresponds to 
/o ~ (J o(R/Ro) 2 ~ k j or to a radius i?i that can be defined by <r M (i?i) = 1 and is given by 

, n \ l/(2-fc) 

i«i ~ #o (^J • (28) 

For k = 2 both T CD = T u and <tcd = do not change with radius, so that generally a u 
is either always below 1 or always above 1, corresponding, respectively, to regimes I and II 
that are discussed below, so that in this case there is no radius R\ where a u (Ri) = 1. 

If the magnetization in region 3 or just behind the CD is low, a u ~ ctcd *C 1, then 
T CD « 2cr according to Eq. (j2T|) . so that Eq. f )25|) implies 

UJ • ™ E ~UfJ UJ <<:1 - (29) 



If, on the other hand, the magnetization in region 3 or just behind the CD is high, <j u ~ 
o"cd 1) then Eq. ( 12"7|) implies r 4 « l/2p (since a > o u in all of the region behind the CD), 
and in particular Tcd = r.4(£«) ~ l/2p u , so that Eq. (|25|) gives 

"""U^J UJ ■ a " KffCDK llAj UJ >>: ' 

1/4 / £ \ ( fe - 2 )/ 4 / p \ (*-2)/4 

V 32a J J ~ Fcr 

where i? cr ~ -Rc^cr is the radius at which T C d reaches the value r cr when a u m o"cd 3> 1, 
and an expression for this radius is provided in Eq. (IIS]) below. 



r m „(i^r(£Y™"~r.(#| ,30) 



4.2. Regime I 

From the derivation above it becomes clear that for / ^> cr 7 Q ~ 2k the external medium 
would hardly affect the acceleration phase, and the magnetized shell would accelerate es- 
sentially as if it were expanding into vacuum (as described in paper I, and summarized in 
§ [3]). This can be seen from the fact that this condition corresponds to a u (R c ) <C 1, i.e. 
that even by the coasting radius R c ~ Rq(Jq the region of the original shell that had been 
affected by the external medium (region 3) occupies only a small part of the flow near its 
head that carries a small fraction of its energy. The transition, where fo ~ &o~ 2k , corre- 
sponds to the equality of the coasting radius (or distance), R c ~ RqO'q, and the deceleration 
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radiuq^j. In planar symmetry with a constant external density p\ (which corresponds to 
k = 2 in spherical symmetry), conservation of energy implies E ~ l a p c 2 ~ (y^ldecpic 2 and 
thus the deceleration distance is given by /dec ~ /oPo/°"oPi = ^o/o/co (where for simplicity 
we discard factors of order unity) so that indeed /dec ~ h ~ ^o " 2 corresponds to f ~ (Tq, as 
it should. For spherical symmetry, energy conservation reads 22 ~ 2?qctoPoc 2 ~ cTqAR^c 2 , 
implying a deceleration radius 22 dcc ~ (E / 'a^Ac 2 ) 1 ^ 3 ^^ ~ 2? (/ /o r o) 1 ^ 3 ~ fc \ so that indeed 
2?dec ~ 2? c ~ 2? CTo corresponds to f ~ cTQ _2fe or A ~ R^p aQ k ~ 7 . Note that in this regime 
2? d ec essentially corresponds to 2? r that is given in Eq. §5§ where T is replaced by T(R C ) ~ cr . 

For k < 2, a u increases with radius (see Eq. |29] ) and since we have seen that regime 
I corresponds to a u (R c ) <C 1 this implies that o u <C 1 all along. For 2 < k < 10/3, on the 
other hand, a u decreases with radius passing through the value of 1 at a radius Ri given 
by Rx/R c ~ (/o/o"o~ 2fc ) 1/(2 ~ fc) < 1- This would be physically interesting only if R x > R , 
which corresponds to crl~ 2k <C fo <C erg. In this parameter regime Ro < R\ < R Cl so that 
Eqs. ([29]) and ([30]) imply that a u « cr CD ~ (R/R^ 2 -^^ > 1 at R < R < Ri while 
o'u ~ o"cd ~ (R/RiY 2 ~ k ^ 2 < 1 at 2? > 2?i (or at 2?i < 2? < 2? cr , as we shall see later). For 
both < 2 and 2 < k < 10/3 we have cr u £s cr CD <C 1 at R > 2? c . 

One can find the time when the reflected rarefaction wave reaches region 3, = £ u , or 
the corresponding radius R u . Relying on the derivations in paper I, one obtains R u /R c ~ 
2c7u 3//4 , which upon substitution of <j u (R u ) from Eq. fT29|) and solving for 2? u gives 



24 2 / 3/ V 4 - 3fc (Rns Xli - 3k 



Rc aT{Ru) \^aol- 2k ) \R 



» 1 , (31) 



where 2?r§ (discussed below) is the radius where a strong reverse shock develops. Therefore, 
clearly R u < Rus, and the reflected rarefaction reaches region 3 well before a strong reverse 
shock develops. The rarefaction wave also reaches the CD (£* = £cd at a radius 22*,cd) 
within a single dynamical time from reaching £ u (i.e. 22*cd ~ R u ), 

A# _ Ai = £ CD - & = ~ i-fep/Wfep) = l + /3cD/3ms(/3cp) x f32) 
R u ~ t u /3*(/3cd)-/3cd feD+ftn B (/3cp) _ ^ i _ /3 CD /3 ms (/3 CD ) ~ ' 1 j 

l+PCPPms(PCP) ' 

1 /2 

where A2? = R^cd—Ru and the last approximate equality is valid since /3 ms (/3cD) ~ cr u <1. 
Once the right-going rarefaction wave reaches the CD, this triggers a gradual deceleration of 
the CD, which is initially weak as the rarefaction is weak at this stage since w ms = Jou <C 1. 



13 The deceleration radius i?dec is the radius at which most of the energy is transferred to the shocked 
swept-up external medium. Here its post-shock Lorentz factor is T 2 ~ Co and therefore the energy given to 
a swept-up external rest mass M is (r| — l)Mc 2 ~ <r 2 Mc 2 , and i?dcc is given by E ~ a\M(r < i?dec)c 2 . 
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In regime I, which corresponds to fo 3> oJ^ 2k -<=>- a <C T CJ - ~ (/oOo) 1 ^ 8-2 ^ or R c <C 
Rdec ~ -Rr> there are three main stages in the dynamics of the shell (see Figs. [2] through 
[7]): (i) initially (at Rn < R < R c ) the shell accelerates, its typical Lorentz factor increasing 
as (r) ~ (a R/ Rq) 1 / 3 while its typical magnetization decreases as (cr) ~ ^^(R/Rq)^ 1 ^ 3 
(since magnetic energy is converted into kinetic energy while the total energy is conserved, 
(r)(cx) ~ o"o), (ii) at the coasting radius, R c ~ RoPq, the kinetic energy becomes comparable 
to the magnetic energy, (cr) ~ 1, so that at R c < R < Rd ec most of the energy is already in 
kinetic form and the shell coasts at (r) ~ while its magnetization decreases as (cr) ~ R c /R, 
(iii) at .Rdec ~ {E / a 2 Ac 2 Y[ {3 -^ ~ i? r (r -> cr ) most of the energy is transfered to the 
shocked external medium! 14 !, and at R > Rdec the flow approaches the BM76 self-similar 
solution where (r) ~ (E / Ac 2 ) 1 ^ 2 R^ k ~ 3 ^ 2 . This is summarized in the following equation: 




Rn < R < R, 



C ) 



R c < R < R, 



dec 



' on(R/R c ) 1 ' 3 R <R<R C , 

(r) ~ { cr R C <R< R dcc , (a) - 

, <J (R/Rdec)~ R > Rdec , 

(33) 

Please note that in regime I, T CD ~ a at R > max(_R > Ri) R c *C i?dec- However, at 
i? > i?dec as the original magnetized shell becomes part of the BM76 self-similar solution 
its Lorentz factor is ~ Tcd and it decreases with time along with its magnetization and 
total energy. As long as it is relativistically hot and thus part of the BM76 solution, its 
Lorentz factor scales as Tcd oc R( 2k ~ 7 ^ 2 oc T^ 2fc_7 ^' 4 ( 4 ~ fc ^ while its magnetization decreases 
as cr oc R( 2k ^ 9 )/ 2 oc 2 1 ( 2fc - 9 )/[ 4 ( 4 ~ fc )] ; where T ~ R/cTq D is the time when radiation from 
the original magnetized shell reaches the observer. However, since the reverse shock is only 
mildly relativistic the shell's temperature quickly becomes sub-relativistic and it deviates 
from the BM76 solution (and the corresponding scalings above), decelerating more slowly 



fjKobavashi & Sari!l2000f ). 



In regime I, the typical magnetic pressure in the ejecta shell at R c is p m (R c ) ~ p(R c )c 2 ~ 
PocVq -5 (where p is its typical or average proper density), while the pressure of the shocked 
external medium is P2(R C ) ~ Pi(R c )c 2 (T 2 = pi(R )c 2 aQ~ 2h , so that the typical or average 
magnetic pressure in the shell is much larger, p m (R c ) / P2(Rc) ~ /o/ cr o~ 2fc ^> 1- However, at 
larger radii the two pressure scale as p m oc R~ 4 and P2 oc R~ k so that their ratio drops with 



14 At R < Rdec ~ Rr the shocked external medium holds only a small fraction of the total energy, 
E cx t/E ~ (i?/i?r) 3 ~ fc < 1 (for k < 3 for which the forward shock decelerates). 
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radius as p m /P2 oc R and the two pressures become comparable at -Rrs, where 

-RRS ~ Rc y a 7-2k J ~ -Rcr , -Rdec ~ R C J " (^ 4 ) 

A strong reverse shock must form at i? ~ -Rrs, since at that stage the magnetic pressure 
can no longer balance the thermal pressure of the shocked external medium at the CD, and 
a new source of pressure is needed, which comes in the form of thermal pressure that is 
generated by the reverse shock that develops and soon becomes dominant. While a weak 
reverse shock might develop earlier, at R < -Rrs the thermal pressure it generates would 
be much smaller than the magnetic pressure, so that it would not have a significant effect 
on the dynamics and would dissipate only a small fraction of the total energy. The reverse 
shock is initially Newtonian, until it becomes mildly relativistic at -Rdec- This can be seen by 
balancing the pressure behind the forward shock, p 2 ~ pi(R)c 2 a 2 ~ pi(R )c 2 &l~ 2k (R/ R c ) k , 
with the (predominantly thermal at R > -Rrs) pressure behind the reverse shock, pus ~ 
p{R)c 2 u\ s ~ p c 2 a ~ 5 (R/R c )~ 3 u 2 ls , which implies a reverse shock upstream to downstream 
relative 4- velocity of u RS ~ (R/R dec ) ( - 3 ^ k ^ 2 - 

This is the familiar "thin shell" case for the deceleration of an unmagnetized initially 
coasting shell (described in §|5J). The shell starts spreading significantly (in the lab frame) at 
R c ~ _Ror 2 (i? c ) ~ RqUq., resulting in the formation of a reverse shock that becomes thermal 
pressure dominated around -Rrs, and gradually strengthens until it becomes mildly relativis- 
tic near its shell crossing radius, which is the deceleration radius, .Rdec ~ (E / a 2 Ac 2 Y^ 3 ^ k \ 
Near -Rdec, where most of the energy is given to the shocked external medium, and where 
the reverse shock crosses most of the shell, the typical magnetization of the shell is low, 
(a) ~ R c /R dec ~ (a /r cr ) 2 ( 4 - fc )/( 3 - fc ) ~ (a 7 - 2 V/o) 1/(3_fe) < 1 (where I have identified A in 
Eq. [TT] with Rq). Note that this regime corresponds to T(R C ) ~ a® T cr , which can also 
be expressed as r cr /a ~ (/oOo* -7 ) 1 ^ 8-2 ^ ~ Y^~ k ^( 8 ~ 2k ^ ^> \ where in the expression for 
T (Eq. [8]) one substitutes A — > R and r — > <J , thus clearly corresponding to the un- 
magnetized (or low magnetization) thin shell case. A larger magnetic field downstream (and 
also somewhat upstream) of the reverse shock is possible due to magnetic field amplification 
in the reverse shock itself, which may allow for a reasonable radiative efficiency coupled to 
the rather effective energy dissipation in the mildly relativistic reverse shock. 



' .Rr, and 

n„\— /?,,. fn,i,A~ f n >\~ f^crY ( fo N 1/(4_fc) 
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Rr J Rr V R, 
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2-h 







7-2k 
°0 



> 1 . (35) 



The ordering of the relevant critical radii in the different regimes are given in Tables |2] and 
|3j In regime I with k < 2 or with 2 < k < 10/3 and fo > o"q we have -Rq < R c < R u ~ 
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-R*,cd < -Rrs ~ Rcr < Rdec ~ -Rr while for 2 < < 10/3 and Ug 2k < fo < <Jq we also have 
the critical radius _Ri so that _R < R\ < R c < R u ~ -R*,cd < -Rrs ~ -Rcr < -Rdec ~ -Rr- 

4.3. Regime II 

This regime corresponds to a^ 3 <C fo <C lTq~ 2A: -<=>- 0q 12 3fe ^ <C r cr <C cr -<=>• T cr <C 
do "C Tc/^ 12 3fe ' ) , where the condition r cr <C cr corresponds to i?dcc ~ -Rcr ~ -RcTcr <C 
Rq<Jq ~ i? c . As we shall see below, this also implies that a,!/ 3 <C T C d(R u ) <C a and 
1 < (7„(i4) < cr 2/3 . 

For k < 2, <j u increases with radius, and since in regime II it is larger than 1 at 
R u , it passes through the value of 1 at a smaller radius Ri that is given by R\/R c ~ 
(,/o/oo~ 2fc ) 1/(2 ~ fc) < 1, and the ordering of the critical radii is R 1 < R u < R CI < R c . As 
in regime I, also here in regime II, R\ is physically interesting only if Ri > Rq, which now 
corresponds to <7q < f < o]^ 2k . In this parameter range a u ~ o"cd increases with radius as 
o u « ct C d ~ (R/Ri) {2 - k)/2 < 1 at Rq < R < R x and as a u « ct CD ~ (R/ Ri) {2 - k)/A > 1 at 
Rx < R < R cr . For cr 1/3 < f < a% we have R x < Rq, and a u « cr CD ~ (R/ RjW* > 

1 all along. Altogether, for k < 2 we have cr u « <tcd ~ (R/RiY 2 ~ k ^ 2 < 1 at i? < 
i? < min(i? ,i?i) and as cr u « ct C d ~ (R/Ri) {2 ~ k)/A > 1 at min(i? ,-Ri) < R < R CT - For 

2 < k < 10/3, on the other hand, R x /R c ~ (iV^cr) (4 ~ fc)/2 (/o/^ 2fc ) V(2 ~ fc) > 1 so that 
R x > R c > R CT and cr u « <x cd ~ (i?/-Ri) (2 ~ fe)/4 > 1 all along. 

For = 2 we have a self-similar solution for the rarefaction wave, thanks to the equiva- 
lence of the cold MHD equations for a spherical flow to those for a planar flow, which make 
it easier to explicitly calculate much of the relevant dynamics. For a general value of k we do 
not have this privilege, and I have relied on the approximation that this self-similar solution 
still approximately holds in this case where Tcd = and £ u gradually change with time. 

In order to further justify this, I now provide an alternative derivation of Eq. fl30l . The 
pressure balance at the CD reads (-E>cD/rcD) 2 /87r m a(4/3)r^ D piC 2 , implying 

^"G&f*™ 4 - (36) 

where Lqd ~ C-B^d-Rcd * s the instantaneous Poynting flux through a static spherical surface 
at r = -Rcd- Note that Lcd is close to the mean (isotropic equivalent) luminosity (or 
power) of the source, L w Ec/2R w irp aoc 3 Rl [identifying the initial width of the shell 
A with its initial radius Rq, where the shell initially occupies the region < r < R , 
while E £s E EMt0 = 2iip a c 2 RQ and p — p(t = 0,r = Rq)], only where the magnetization 
parameter just before the CD is large, cxcd 3> 1, which corresponds to Tcd "C er . In this 
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case Lqd ~ £ and we have 



CD 



V32avrAc 3 J ~ \Q4a7rAc 2 R ) ~ V 32a / 



R 



R \ (fc-2)/4 / \ ( fe ~ 2 )/ 4 / # \ ( fc " 2 )/ 4 



r S- l,/2 (£J , (37) 



2rcD V 3/o / \Ro ) \Ri. 

This is valid as long as the value of the lab-frame magnetic field B at the CD (i.e. the 
head of the outflow) is close to its original value, i.e. for <tcd = °"« ^ 1, which holds at 

max(i? ; R\) < R < Rdec- 

The condition that fo ^> a^ 3 in regime II implies that Tcd(-Ro) ~ Icr - *^ 2 ~ (jo^o) 1 ^ 4 3> 
er , and therefore at to region 4 (see Fig. [1]) holds most of the volume and energy, and (r) ~ 
Gq Z (R/ Rq) 1 ^ ~ ((Tot/to) 1 ^ 3 at t > t - At this stage the typical or mean value (weighted 
average over the energy in the lab frame) of the Lorentz factor within the magnetized shell, 
(r), increases with time, while for k = 2 the Lorentz factor of the uniform region at its 
front, r(£ u ) = T 3 = Tcd, remains constant. More generally, Tcd(-Ro < R < -Rcr) is given by 
the smaller between the expression in Eq. fl37|) and 2o"o. This acceleration (increase in (T)) 
lasts until the secondary (or reflected) rarefaction wave finishes crossing region 4, i.e. until 
£, « 1 - 2i? /ct « 1 - 2^/^ equals 



1 /o \ 3/4 3(2-fc) 

( (3 C B - /3 ms (/3 CD ) ^ 1 - (r ms /r CD ) 2 ^ a _ 1 /'32aa 1/3 \ { R\ * 



1 - /3 C D/3 ms (/3 C D) i + (r ms /r CD ) 2 r* D i 3/0 7 V#o 



(39) 

(where since r ms ^> 1, we have r 2 ^ « u 2 ^ = a = eroP ~ c"o/2r and specifically T^^/^cd) ~ 
ct /2Tcd), at R = R u , which corresponds to 



Ru I 3/ \ /-nWnN f (T oRu\ ( 3/o<7, 



1/3 /op 3-fc 



<r> W ~ ( ^ ) « ( ) « f cd(^) = r 



10-3fc 



This implies 



(40) 

R c ~ \2^aal 2k ) ~ 4a*Mj ^ 1 ' (41) 
which is different from the result for regime I (see Eq. [2D), where R u /R c « 2a u 3/4 (R u ) > 1. 



-Ru / 3/o 



10-3fc 



-23 - 



At this stage (R = R u or t = t u ) most of the energy in the flow is in_| region 3, which 
moves with T 3 w T CD given by Eq. f )37|) . that represents (T) at this stage (R u < R < -R*,cd)- 
Region 3 is gradually crossed by the right going rarefaction wave, until it reaches the CD at 
R*,cd ~ Rcr ~ -Rdec (as shown in detail below), which marks the end of this stage. At that 
point most of the energy is in the shocked external mediuu@, and the flow approaches the 
BM76 self-similar solution (first the rarefaction wave crosses region 2, within a few dynamical 
times0 and then the adiabatic BM76 self-similar solution is quickly approached). 

The width of region 3 at t u (when the rarefaction wave reaches £ u ) in the lab frame is 
A 3 = ci M (£ C D - £«) ~ R u ^cd( R u) ~ 2-Ro- In region 3, 

n _ ficD + /3ms (Pcd) _ j 1 ~ 1 _ 1 C42) 

i + /3cD/3 ms (/3c D ) ~ 8r 2 CD rL(/3 CD ) ~ 4a r CD ' 1 ' 

so that (1 - /?*) < (1 - /3 C d) ~ V 2r cD and therefore Aw = (/3„ - /3 C d)c « c/2r^ D and the 
increase in radius, AR = R* t cD ~ Ru (or time, At = t^cD — t u ~ Ai?/c), during the time it 
takes the rarefaction wave to cross region 3 is Ai? 2Rqc/Av ~ 4i?orcD ~ -^dec for k = 2, 
while more generally 



**,cd rR*,cn ^jj> jj> rR*,CT>/Ru 



2— fe 



2 



(4-fc)/2 



4fl ° r g D(fi " ) « « tocffl.) » 1 



15 This can be seen as follows for k — 2. The pressure is continuous across the CD, and therefore the energy 
density of regions 2 and 3 in the lab frame is similar, and their relative energy is determined by their relative 
width in the lab frame. For region 2, using the uniform velocity approximation £ s h — £cd = /? s h — /3cd ~ 
l/2r2 D - l/2r s 2 h a l/4r2 D (for the BM76 solution £ sh - £ CD « (1 - X CD )/ 2r cD « 1/4-601^, which is 
rather similar), while for region 3, £cd — £« ~ a o/^cD ~ V^cd — (^o/^cd — 1)/2Fq , and therefore the 
width of region 3, A3, is larger than that of region 2, A2, by a factor of A3/A2 ~ 2(2tro/rcD — 1) 3> 1, since 
Tcd *C (Jq in this regime. 

16 At R u < R < i?dcc ~ -Rcr only a small fraction of the total energy is in the shocked external medium, 

E c ^ t /E~{R/R c ^- k )/\ 

17 For k = 2, making the approximation that the region between the CD and shock front has the constant 
velocity of the CD and that r s h = V2Tcd and using Eq. ([2T]) one obtains that during the time the rarefaction 
wave travels from the CD to the shock front the radius increases by a factor of 1 + (1 - Xcd0 2 (v / 3 + 1) /(3 - 
•\/3) ~ 2.87. If we self-consistently use the above assumption to estimate the width of this region (even 
though this is not fully self-consistent) this gives (R s h — Rcd ) / Rcu ~ 1/ 4Fcd instead of Eq. (|2"TT) , and a 
growth in radius during the rarefaction crossing by a factor of 1 + ( yo + 1 ) / (3 — v3) ~ 3.15. In both cases it 
is close to a factor of ~ 3. This factor is relatively large since the sound speed in region 2 in "only" c s c/ yo 
(as it is unmagnetized but relativistically hot, while regions 3 and 4 are cold but highly magnetized) and the 
shock front moves somewhat faster than the fluid in region 2. 
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[2(4-k)a CD (R u )} 



2/(4-fc) 
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2/(12-3fc) 
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Ru 



(43) 



so that the rarefaction reaches the CD at .R* cd ~ A.R ~ R CT ~ -Rdec- 



The deceleration radius in this regime can be obtained by equating the initial mag- 
netic energy to the energy of the swept-up external medium, E « 2n p Q a c 2 R^ ~ [47r/(3 — 
k)}Ac 2 R 3 d ; c k T C B{Rdec) 2 , which implies 



R 



dec 



-Rn 



(3 - A;) 2 8a/ cT 



-, 1/(4-*) 



~ r CD (i? dcc ) ~ 



r 



R cv 
R 



and therefore i?, 



dec 



i? cr where _R cr is the radius at which Tcd = Tbm, 
'(3- k)E^ 1/2 



and 



BM 



Ait Ac 2 



R -(3-k)/2 



(44) 



(45) 



is the typical Lorentz factor during the subsequent constant energy self-similar (BM76) stage. 
Estimating the value of R CT from Eqs. ( 13?)) and f|4"5]) and identifying _R with A gives, 

1 i /( 4- fe ) f 9.3 x 10 16 aV4 C -i/4 no 1/4^1/4^/4 cm (fc = Q) ^ 

(46) 

5.3 x 10 15 a 1 /2 C -i/2 A -V2 jE; i/2 T V2 cm (jfc = 2) s 



4(3 - k) 2 aEA 
3ttAc 2 



or 



Rr 



1 fER 



Ro Rq V ^c 2 



l/(4-fc) 



(/o^o) 



1/(4-*) 



(47) 



During the initial acceleration (at -R > i?o), (r) ~ (a R/ Ro) 1 ^ 3 ■ This lasts until most 
of the energy is transfered to the part of the magnetized shell with T ~ Tcd, which occurs 
at a radius R u , Lorentz factor T u , and magnetization a u (R u ) given by 



Ru f, _ 1/3 \ 3 /( 10 - 3fe ) 



r ~ 



<T0 



Cu(Ru) 



//3 



~ I /O^O 17 



/3 x l/(10-3fc) 



Moreover, 



r« ~ r cr ( — 

O"0 



(fc-2)/(10-3fc) 



r 



cr 
O"0 



2(4-fc)/(10-3fc) 



(4f 



(49) 



so that near the transition to regime I, T u ~ r cr ~ cr and R u ~ _R cr ~ i? c 



In Regime II we have 1 <C /qCq 1/3 ^ a o^ 10 3k ^ 3 , which corresponds to 1 R u / R <C cr 2 . 



(i.e. Rq <C -R u *C -R c ), o"o /3 < T u < (T and 1 <C a u (R u ) <C (To 7,3 . The different critical radii 
are related by 



2/3 
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so that for 2 < k < 10/3 or for k < 2 and a^ 3 < / < o"q we have R\ < Rq < R u < 
R cr ~ Rdcc ~ R*,cd *^ -Re &nd -Ri is irrelevant (as cr u 3> 1 all along), while for k < 2 and 
<7q < /o < o"o _2/c we nave Rq < Ri < R u < R C t ~ -Rdcc ~ -R*,cd < R c an d -Ri is relevant. In 
all cases R c is not relevant physically (since it looses its meaning as a coasting radius). 



The typical magnetization of the shell in the intermediate stage is 

(<T)(Ru < R < R CT ) ~ <7CD 



cro ( R\ {2 - k),A cro f R^ {2 - k)/A 



Tcd r« V-R« 



Tcr V -Rcr 



(51) 



so that at i?dcc ~ Rcr we have (cr) ~ o"o/r cr 1- Thus, altogether in regimes I and II we 
have 

(cr /r cr ) 2 ( 4 - fc )/( 3 - fc ) ctq < T cr (regime I) , 

(52) 

o"o/r c r o"o > T cr (regime II) , 



(a) (Rdcc) ~ 
while in regime II we have 



{aoR/Ro) 1 / 5 ~ T^R/Rtfl* R < R < R u , 



(T)(R) ~ ^ 



r cr (R/R CT )( k -w 



R u < R < R cr , 



-R > -Rcr , 



(53) 



cr Rq < R < max(i? 0) -Ri) , 

r CD ( J R) ~ < (54) 

aoiR/R^-W 4 ~ r cr (/?//?„) (fe - 2)/4 max( J R , R^ < R < R cr , 



a 2/3 (R/R 



-1/3 



(cr /r u )(i?/ J R u )-V3 R <R<R U , 



(a)(R) ~ < (55) 
(cr /r u )( J R/ J R tt )( 2 - fe )/ 4 ~ (cr /r cr )( J R/ J R cr )( 2 - fe )/ 4 R> u < i2 < R> c 



{R/R^-W ~ (cT V/ ) V2 (jR/jRo) (2^)/2 ^ < < max(jR()) ^) s 

cr CD (i2) ~ <( (56) 
(R/R!)*- 2 -^ 4 ~ (cr /r cr )(i?/ J R cr )( 2 - fe )/ 4 max(i? , R ± ) < R < R cr . 



4.4. Regime III 

In regime II we had r CD (_R ) 3> °"c/ 3 so that the plasma near the CD was super-fast- 
magnetosonic with respect to the "wall" already at i? ~ i? , with r CD /r ms (2rc D /cr ) 1 / 2 3> 
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1, and thus not in causal contact with the source. Here, in regime III, we consider what 

1 /3 

happens when 1 <C r CD (_R ) <C a ' . In all the regions behind the CD, the fast magnetosonic 
Lorentz factor is given by r^ s w = a > cr CD « cr /2r CD 3> 1, so that as long as 
Tcd *C c^ 3 the flow remains in causal contact with the "wall" or central source, r/r ms < 
(2T 3 CD /a ) 1 / 2 < 1. Thus, the flow remains roughly uniform and the conditions just behind 
the CD are representative of the typical values in the shell, (T)(R ) ~ Fcd(-Ro), implying 



(T)(Ro) - It-V 2 ~ al'* (W*) = (/o-o) 1/4 , (57) 
(a)(R ) ~ a rr 4)/2 ~ °T (fo^y 114 = (k^V 11 ' , (58) 



so that the conditions near R cr are very similar to those in regime II: Rd ec ~ -Rcr ~ -Roller , 
(T)(R CV ) ~ T cr , (a)(R cr ) ~ (a)(i? )( J R cr /i? ) (2 - fc)/4 ~ <T Q /T a , and r CD (i? dec ) ~ (r)(i? dcc ) ~ 

1 /3 

r cr . This implies that regime III, defined above through the condition 1 -C Tcd(-Ro) ^ 0"o > 
corresponds to a x / ^ Co > 1 ^ ^cr < a- 2 /( 12_3fc ) or ct > p(i 2 ^ 3fc )/ 2 ^ \ 

In this regime Tcd is still given by Eq. (|36|) while 

t /3cd - /3 ms (ficp) ^ 1 - (r C D/r ms ) 2 ^ 1 ^ 4T|p _ _ 

1 - /3 CD /3 m s(/3cD) ~ 1 + (r CD /r ms ) 2 ~ ! + H£|d ~ a ~ ' 1 J 

so that region 3 initially occupies most of the volume, A 3 (t < t u )/ct = £ CD — £ u ps 2 — 
l/2rQ D + 4r 3 JD /o" ~ 2. This demonstrates again that already at t = t the conditions 
just behind the CD (region 3) dominate the average values over the original magnetized 
shell, so that (r) ~ T CD and (a) ~ ctcd are given by Eqs. (13"7|) and (I3"8"j) . respectively. 
Region 4 occupies only a small fraction of the total volume already at t — to, A4(to)/cto = 
£« + /3ms,o ~ 41^/00 <C 1, and it is very quickly crossed by the right-going rarefaction wave, 
which reaches region 3 (£* = £ u ) at i = t u (and -R = _R U ) that corresponds to 



hi _ 1 ~ ^fa) w A 4fa) w 2r cp < x _ (qq) 
to \fi*(to) - £u]ct 2ct cr 



This implies that R u ~ Rq and the time since to when the right-going rarefaction wave 
reaches the CD (i CD — t ) is dominated by its propagation time through region 3. For k = 2, 

^cp_ _ A 3 (t ) w 4 r2 ^ 1 rgl N 

to [/3*(/3cD)-/3 C D]ct 4icD>>1 ' 

so that it reaches the CD at a radius R*,cd ~ c ^cd ~ 4rQ D i? ~ R CT ~ -Rdec- Similarly, 
since in this case R u ~ i? then for a general k- value Eq. (j4"3"j) implies -R*,cd/-Ro ~ AR/R ps 
[2(4 - A;)r 2 D ( J R )] 2 /( 4 " fc ) ~ T 2 r ~ iW^o, so that again, ^, CD ~ Rcr ~ Rdec 
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The effect of the external medium in this regime is very large, in the sense that it causes 
most of the energy to be in the uniform region 3, with a sub-fast-magnetosonic speed relative 
to the "wall". Nevertheless, since this region is still relativistic, it takes the rarefaction wave 
that is reflected from the wall a long time to cross this region in the lab (or wall) frame, and 
this occurs at a large distance from the wall, -R*,cd ~ -Rdcc ~ -Rcr ~ -Rol^r ^ -Ro; near the 
deceleration radius where most of the energy is transferred to the shocked external medium. 
Altogether, in regime III we have R ~ R u < R CT ~ Rd cc ~ R*,cd < R c (see Table |2]), and 
similar regime II, here as well R c does not have a physical significance (and the same also 
holds for Ri, since we always have o u ^> 1 in regime III). 

In regime III, (T)(R > Ro) becomes independent of ao while (c)(R > Ro) scales linearly 
with cr , when fixing L, A, k and i? (which fixes r cr ) while letting a and p ^ l/°o var Y 
(since a ^> 1 we have L ~ Ec/Rq ~ (JqPqC'Rq oc <7 p , so that fixing L implies that 
Po oc l/o"o). Such a variation of the parameters means fixing the overall properties of the 
flow and changing only its composition or magnetization (as is done in Figs. |2] through [7]) . 
In this regime the global dynamics become insensitive to the exact composition. This can 
be thought of as the high magnetization limit, where the behavior of the outflow approaches 
that of an electromagnetic wave that is emitted at the source and reflected by the CD, where 
the time when the back end of the finite wave reflects off the CD corresponds to the time 
when the right-going rarefaction wave reaches the CD, ten ~ -R*,cd/c 

Alternatively, as is done in Figs. Eland El one could fix the properties of the magnetized 
flow: L, R , o" , po (and thus also R c ~ Ro^o) anc ^ var y the normalization of the external 
density: A or pi(Ro) = AR^ k (while fixing its power-law index, k), which effectively varies 
r cr and R cr . It can be seen from Figs. [8]and[9]that as the external density goes to zero we have 
fo — > oo, r cr ~ (/o (T o) 1 ^ 8 ~ 2A '' ) — >■ °° an d -Rdec ~ Rr — > oo, and this the solution approaches 
that of expansion into vacuum (or the extreme limit of regime I). As the external density 
increases f , r cr and i?d ec all decrease, until when r cr ~ a , f ~ o"o~ 2fc an d -Rdec ~ -Rr ~ -Rcr 
there is a transition to regime II. As the external density increases even further a transition 

2/(12 3k) 1/3 

to regime III occurs when r cr ~ a , fo ~ cr and R u ~ Rq. Finally, when the 

external density becomes so large that r cr ~ 1, f ~ a^ 1 , and R CT ~ -R 0) the flow remains 
Newtonian and there is a transition to regime IV that is discussed below. 



4.5. Regime IV 

For a sufficiently large external density, fo ctq -1 , the formal expression for r cr gives 
r cr 1 and the flow remains Newtonian. If we consider a source that is active over a time 
to then when the central source finishes ejecting the highly magnetized outflow, it would be 
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bounded within -Rcd(^o) ~ PcT>{to)cko ~ /9cd(^o)-Ro, where I neglect factors of order unity 
for simplicity (here Rq is still defined through the relation R w ct Q , even though it loses 
its physical meaning from the relativistic regime). More generally, at t < t the radius of 
the CD satisfies Rcu{t) ~ /3cd(^)c£ ~ f3cD(t)Rot/t . For a tangential magnetic field, which 
scales as B/B Q ps Rq/R, the magnetic pressure at -Rcd, ~ -S 2 [_Rcd(^)] ~ [Ro / Rcvif)] 2 B 2 ~ 
<?oPoc 2 {t/to)~ 2 (3c^(t), would be balanced by the ram pressure of the shocked external medium 
at the frame of the CD, ~ p 1 [ J R CD (t)]/3^ D (t)c 2 ~ p 1 (R )c 2 (t/t )- k /3l~ k (t), leading to 



This implies /3cd(^o) ~ (/o^o) 1 ^ 4- ^ ~ T c 2 r <C 1, which demonstrates self-consistency by 
showing that the flow is indeed Newtonian in this regime. The magnetic energy in the 
original outflow, at R < R CD (t), is given by E B [R < Rcn(t)] ~ Rc D (t)B 2 [R CD (t)] ~ 
(t/t )(3 C D(t)R 3 B 2 ~ {t/t )pcn{t)E « C D{t)L o t, where L « Sgi^c and £ « £o*o are 
the injected luminosity and corresponding energy over a time to for a relativistic out- 
flow. This would violate conservation of energy, if the outflow emanating from the cen- 
tral source was indeed relativistic (this is basically the well-known <7-problem)0 How- 
ever, since in this regime the outflow is sub-sonic (or sub-fast-magnetosonic) and Newto- 
nian, the information about the existence of the external medium must propagate back to 
the source producing a back reaction that results in a Newtonian outflow with a speed 
~ (3cr>(t). For such a Newtonian magnetized outflow the electromagnetic luminosity is 
L = 47tR 2 (c/Att)\E xB\= B 2 R 2 (3c } and for (3(t) ~ /3 CD (t) this gives L(t) ~ p CI) {t)B 2 R 2 c « 
f3cjy{t)LQ and E(t) ~ tL(t) ~ (3cD(t)L t, which is consistent with the above estimate. Even 
during the initial injection phase (t < t ) the shocked external medium holds a good fraction 
of the total energy at any given time. After the injection stops, at t > to, most of the 
energy is quickly transfered to the shocked external medium on the dynamical time (up to 
t ~ 2t or so), and the flow settles into an adiabatic Sedov- Taylor solution with velocity 
P(t > t )c ~ [£7/A* 3 -*] 1/(5_fe) ~ /Wto)c(t/t o r (3 " fc)/(5 - fc) , radius R ~ (Et 2 / 'Afl^ and 
energy E = E(t ) ~ /3 C d(^o)^o- 

If we start with a magnetized spherical shell or "ball" of radius Ro, initially at rest, then 
in this case the magnetic pressure at Rqd ~ R® is ~ B 2 ~ UoPoC 2 and equating it to the ram 
pressure of the shocked external medium, ~ p 1 (i? )/3 2 , D c 2 , implies /?cd ~ (/o^o) 1 ^ 2 ~ r^7 fc <C 



18 For a steady central source that ejects a magnetized outflow over a finite time t$ this regime resorts 
back to the well known a problem, where in ideal MHD the stored magnetic energy grows quadratically 
with the injection time of the central source, while the actual injected energy grows only linearly with this 
time, implying a breakdown of one or more of the underlying assumptions. This may be relevant, e.g., for a 
millisecond magnetar born inside a collapsing massive star, as a possible progenitor of long duration GRBs. 




(62) 
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1. Therefore, it would significantly increase its radius and transfer most of its energy to the 
shocked external medium on its dynamical time, which corresponds to a timescale of 

tdec ~ -5^- « £~ ~ (/0^0)~ 1/2 t ~ r k c -% » to • (63) 

Pcdc Pcd 

At t > tdec most of the energy is in the shocked external medium and the flow settles 
into an adiabatic Sedov- Taylor solution with velocity f3(t > tdcc)c ~ [E/At 3 ^ k ] 1 ^ 5 ^ k ^ ~ 
/3cD(tdcc)c(t/t )" (3 ~ 5)/(5 ~ fe) , radius R ~ (St 2 /^) 1/(S-fc) and ener gY E = E o ~ ^opoc 2 . 



4.6. Regime II* (10/3 < fc < 4): 

When the external density drops very sharply with radius, k > 10/3, then Tcd initially 
grows with radius faster than R 1 ^ 3 , which has interesting implications. Regime II that exists 
for k < 10/3 disappears for k = 10/3 and reappears for 10/3 < k < 4 in a different form 
that we shall call regime II*, which corresponds to (?J ) ~ 2k <C fo <C cr^ 3 . For 10/3 < k < 4, 
Regime I holds for f 3> a^ 3 and regime III holds for a^ 1 <C fo <C o-J^ 2k . 

In regime I, for / > o"q •<=>- oo < Tcr fe ' )/ ' 2 there in no R\ (as ctcd <C 1 and Tcd ~ o"o all 
along) and the shocked external medium decouples from the magnetized shell at a decoupling 
radius of i?d cp ~ Rq < R c , as the forward shock accelerates down the steep external density 
gradient, and the shocked external medium carries only a very small fraction of the total 
energy, E ext /E ~ er //o < <j q 2 <C 1. In the parameter range a^ 3 < f < (7q •<=>■ rlr - ^ 2 < 
cr < TcJ 2 3fc ^ 2 ; on the other hand, there is a radius Ri where ctcd ~ <5u = 1 and the ordering 
of the critical radii is Rq < R\ ~ -Rdcp < R c , so that Tcd ~ 0"o/°"cd i?( fc ~ 2 )/ 4 following 
Eq. (ETJ until T C d ~ °"o an d ctcd ~ 1 at -R ~ i?i ~ (/o°o" 3 ) > and the decoupling of the 
forward shock from the magnetized shell occurs at -Rdcp ~ Ri- in both cases at R > Rdcp 
the shell accelerates in the wake of the accelerating forward shock, almost as if into vacuum. 
At R > R c it starts coasting and spreading radially, where it can in principle keep coasting 
indefinitely (or more realistically until the assumption of a very steep external density profile 
breaks down, and enough external mass is swept-up that could decelerate the forward shock 
and, in turn, also the shell). 

In regime II* the ordering of the critical radii is Ro < R u < Ri ~ -Rd cp < R c < -R cr 
(see Table |3]). Initially, at Rq < R < R u , the typical Lorentz factor and magnetization 
are similar to those at the CD and are determined by the pressure balance at the CD, 
(r) ~ r CD ~ 00/(0-) ~ o-o/ocd oc i?( fc - 2 )/ 4 following Eqs. §7} and (EE]). At R u < R < R CI 
the bulk of the shell decouples from the CD and accelerates as (r) ~ a / (a) ~ (<JqR/ Ro) 1 ^ 3 
until reaching the coasting radius R c ~ Rqc 2 (where it starts to coast and spread, as in 
regime I), while r CD ~ 0"o/°cd oc i?( fc ~ 2 )/ 4 keeps following Eq. fl5Tj) until T CD ~ a and 
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c cd ~ 1 at R ~ R\ ~ (ioo"^ 3 ) 1 ^ 2 ^- At i? > i?i ~ i? dcp we have ctcd < 1 and r C D ~ o"o, 
while the shocked external medium decouples from the shell, carrying with it only a small 
fraction of the total energy, E cxt /E ~ R\/R c ~ (/o " 2 ^ 7 ) 1 ^ 2 ^ <C 1, as it keeps accelerating 
down the steep external density gradient, with T B m oc i?( fc ~ 3 )/ 2 (where the same energy is 
given to a decreasing amount of newly swept-up external rest mass). The original shell 
keeps coasting and spreading radially in the evacuated region in the wake of the accelerating 
forward shock, essentially as if into vacuum. Near the transition to regime III we have 
fo ~ a o~ 2k an d therefore E ext /E ~ 1, so that the energy in the shocked external medium 
becomes comparable to the total energy. 

In regime III the ordering of the critical radii is Rq < R cr ~ -Rdcp ~ -Rdec, and there are 
no R u or Ri. At R < R < R cr the shell accelerates as (T) ~ T CD ~ cr /acD ~ cr o/(°') oc 
^( fc - 2 )/ 4 following Eq. ( loTj) . where the typical values of T and a are close to those at the CD. 
The shocked external medium decouples from the original magnetized shell at i?dcp ~ -Rcr, 
but since in this case crcD(-Rcr) ~ cro/T C r ^> 1 the rarefaction wave is still strong when it 
reaches R cr so that it effectively decelerates the shell and very little energy remains in the 
original shell at later times (or larger radii), while most of the energy is transfered to the 
shocked external medium, which approaches the BM76 self-similar solution. 

4.7. Summary 

The different regimes are summarized in Tables [2] and [3] as well as in Figure [TO]. Tables [2] 
and [3] provide the ordering of the various critical radii in the different regimes, along with the 
parameter range occupied by each regime, in terms of the initial shell to external medium 
density ratio f = p /Pi(-Ro), initial magnetization a , and T cr ~ (/o^o) 1 ^ 8-2 ^- Figure [TD] 
shows the regions of parameter space occupied by each regime within the two dimensional 
planes spanned by fo-<?o (top left panel for k < 10/3 and bottom right panel for 10/3 < k < 
4), r cr -cr (top right panel for k < 10/3), /o~r cr plane (bottom left panel for k < 10/3). 

In order to gain some intuition for these results and better understand them, it is useful 
to follow the behaviour of the system when varying one key physical parameter and leaving 
the others fixed. First, I vary the initial magnetization ctq (and po oc l/o"o for consistency) 
while keeping fixed the energy (E ~ Lt ps LR /cot L), initial time or length scale (t ps Rq/c 
or R ) and external density (k and A or pi(Ro) = ARq 1 *), and thus also r cr ~ (fo&o) 1 ^ 8 2 ^ 
and R cr ~ RqY^. This corresponds to a constant r cr and varying o~q oc 1/ /o, or in Figure 1101 
to a horizontal line in the top right panel, a vertical line in the bottom left panel, and a 
diagonal line parallel to the fo = Cq 1 line separating regimes III and IV in the remaining 
two panels (showing the fo - o~o plane) . The behaviour in this case is summarized in Figures. [2] 
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through 

In regime I (1 < ctq < Y cr ) the acceleration is almost as if into vacuum: (r) oc R 1 ^ 3 
and (a) oc R~ l l 3 until most of the energy is converted to kinetic form at the coasting radius 
R c ~ Ro<7q, where (r) ~ cr and (a) ~ 1. Then the shell starts coasting (at (r) ~ cr ) 
and its width in the lab frame starts growing linearly with radius resulting in a fast drop 
in its magnetization, (a) ~ R c /R. At R > R c this regime reverts back to the well-studied 
unmagnetized "thin shell" case, with a reasonable spread in its Lorentz factor, 5T ~ (r). A 
reverse shock develops and is initially Newtonian, but strengthens as the shell widens, until 
it becomes mildly relativistic when it finishes crossing the shell at the deceleration radius, 
Rdec ~ Rr ~ (E / (JqAc 2 ) 1 ^ 3 '^ , where the magnetization is low, (cr) (-Rdec) ~ R c / Rd ec ~ 

(t7 /r cr ) 2 ( 4 - fe )/( 3 - fe ) « i. 

In regime II (1 < r cr < cr < r c r 2 ~ 3fc ^ 2 ) the initial acceleration of (r) oc R 1 ^ 3 is limited 
by the external medium at R u , where most of the energy is still in magnetic form ((a)(R u ) ~ 
c"cd(-R«) 3> 1), and thus there is no coasting phase. Between R u and R cr the typical Lorentz 
and magnetization of the shell are similar to those just behind the CD and determined by 
the pressure balance at the CD, (r) ~ &o/(g} oc R^ 2 ^ 4 . A rarefaction wave gradually 
crosses the shell from its back to its front, until reaching the CD at -R*,cd ~ -Rcr ~ -Rdec- 
At that point the shocked external medium starts dominating the total energy and the flow 
approaches the BM76 self-similar solution. 

In regime III (a > r(J 2 ~ 3/c )/ 2 > ;Q the external density is large enough that there 
is no impulsive acceleration stage with (r) oc R 1 ^ 3 . Instead, the pressure balance at the 
CD determines the shell's typical Lorentz factor and magnetization from the very start, 
(r) ~ cro/(c) oc i?( fc-2 )/ 4 at Ro < R < R CT ~ -Rdec, and the dynamics become insensitive to 
the exact composition [i.e. to the value of <7o, when fixing the external density (A and k) 
and the shell luminosity (L) and initial width (-Ro)]- This is the high-o" limit where the flow 
behaves like an electromagnetic wave that is emitted at the source and reflected at the CD. 

There are also more "exotic" regimes, such as regime IV where the external density is 
so high that flow remains Newtonian, or regime II* that exists only for a highly stratified 
external density (10/3 < k < 4) where the external shock accelerates down the steep external 
density gradient and decouples from the original shell, carrying a small fraction of the total 
energy, while the original shell travels in its wake essentially as if into vacuum, similar to 
regime I. Note that regime IV corresponds to r cr < 1 and thus cannot be reached when 
fixing r cr to a value larger than 1 and varying <r oc 1/po- 

A slightly different way to gain perspective about these results is by varying the external 
density normalization (A or pi(Ro) = ARq* 1 or f = p /pi(-Ro)) while keeping the other 
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parameters fixed (k, ao, po, Ro, and therefore also E and L). In this case R c ~ Rqctq 
remains constant while r cr ~ (/oOo) 1 ^ 8-2 *^ an d R C r ~ -^0^ vary, where both / and r cr ~ 
(/o^o) 1 ^ 8-2 *^ f^ 8 ~ 2k ' decrease when the external density increases. In Figure (TU] this 
corresponds to a vertical line in all but the bottom left panel, where it corresponds to a slightly 
diagonal line parallel to the fo = T^~ 2k line (the left boundary of the colored regions). The 
behavior in this case is illustrated in Figures [8] and |9j For a sufficiently low external density, 
corresponding to fo > o-J ) ~ 2k or r cr > o"o we are in regime I, where the expansion is initially 
essentially as if into vacuum, reaching the coasting radius at R c that is independent of fo and 
decelerating significantly only at i?d ec ~ Rr ~ Rc(foo~ 2k 7 ) 1 ^ 3 ~ k ^ oc f^ 3 k \ As the external 

1 /3 

density increases, r CD (i? ) decreases, bringing about first regime II (<r < r CD (i? ) < °o 
or a - 2 ( 3 - fe )/( 4 - fe ) < RJR C < 1), and at even larger external densities regime III (1 < r cr < 
^.2/(12-3*;) QT ^ < h ct /ji < cr^ 4- ^). For the highest external densities (i? cr < -R 0) r cr < 1 
or /o < o-q 1 ) the flow remains Newtonian all along (regime IV). 



5. Comparison with Previous Works 



Sari & Piran 


1995: 


Sari 


1997 



Kobayashi fe Saril l2000l : iKobayashi &: Zhang! 120031 ; iNakar &: Piran! 12004 ) . The main results 



have been summarized in § [2] and extended to a general power-law of the external density pro- 
file, and are consistent with the previous results. The deceleration of a magnetized relativis 



tic shell by an unmagnetized external medium has also been studied (IZhang fe Kobavashi 



2005 



Levinson 



2010 



Lvutikovll201lh . 



Gia nnios, Mimica fe Alovl 120081 ; iMimica. Giannios fe Aloyl |2009| ; iMizuno et al. 



2009 



Zhang fe Kobayashil (120051 . hereafter ZK05) have both considered arbitrary "initial" val- 
ues for the shell Lorentz factor and magnetization, and have attached too much importance 
to the crossing of the shell by the reverse shock, while for a ^> 1 even if a reverse shock 
exists its effect on the global dynamics of the system is very small (it dissipates only a small 
fraction of the total energy, of the order of ~ l/er, and by its shell crossing time only a 
similarly small fraction of the total energy is transfered to the shocked external medium). 
Therefore, the conclusions of that paper are very different from my results. 



Giannios. Mimica fe Aloyl (120081 ) have considered a similar initial setting and argued 



for a different condition for the formation of a reverse shock! 19 !. While the condition for the 



19 Their argument that the shell can be crossed by a fast-magnetosonic wave (and thus come into causal 
contact) faster than by a fast-magnetosonic shock (both starting at the CD) appears to contradict the basic 
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formation of a reverse shock in the ideal Riemann problem addressed in ZK05 is correct, such 
initial conditions are not realistic and the formation of a reverse shock and its properties 
can be sensitive to the exact initial conditions or to fluctuations in the external density, etc. 
Moreover, in the high-a limit even if such a shock exists it has a very small effect on the 
global dynamics, which are the main focus of the present w ork, and therefore this is not 
addressed here in detail. In iMimica. Giannios fc Aloyl ( 120091 ) the problem is addressed with 
a similar initial setup but using high resolution ID R MHD numerical simulation s. The re, 
the regime that is argued to have no reverse shock in iGiannios. Mimica fc Aloyl (120081 ) is 
correctly found to have either a weak or no reverse shock. They also demonstrate numerically 
that the flow quickly approaches the BM76 self-similar solution after the deceleration radius. 



Mizuno et al.l (120091 ) point out that for the Riemann problem of a magnetized shell 



moving relativistically relative to an unmagnetized region (or "external medium") at rest, 
above some critical value of magnetization parameter a there is a rarefaction wave that prop- 
agates into the magnetized shell and accelerates it, and only below that critical value there is 
a (reverse) shock that decelerates the shell. While this observation is correct, this Riemann 
problem is not a realistic setup for the deceleration of magne tized GRB eject a, since it uses 
arbitrary "initial" conditions near the deceleration radius. iLyutikovl (120111 ) has analyzed 
the similar problem of the deceleration of a shell with arbitrary initial Lorentz factor and 
magnetization, concluding that the differences between the magnetized and unmagnetized 
cases are rather small, and involve mainly the existence or strength of the reverse shock at 
early times (which may be non-existent or weak for high magnetizations), rather than the 
global gross properties of the flow. I find that this is a right answer for the wrong question, 
in the sense that the initial setup is too arbitrary to realistically apply to GRB outflows. The 
impulsive acceleration process determines the conditions near the deceleration radius, which 
are therefore not arbitrary, and some regions of parameter space and their corresponding 
dynamical regimes cannot be realized under realistic circumstances. 

Paper I has addressed mainly the impulsive acceleration into vacuum of a highly magne- 
tized shell, starting at rest. However, at the end of its subsection 5.2 it also briefly addressed 
the expansion of such a shell into an unmagnetized external medium. There it has outlined 
the two main dynam ical regimes, which in the current work are referred to as regimes I and 
II. iLevinsonl (120101 ) has also considered the acceleration and of an impulsive magnetized 



notion that a shock must always travel faster than the relevant corresponding wave. It arises since they use 
the formula for the radius at which the reverse shock crosses the shell from Eq. (38) of ZK05 that is valid 
only for a strong reverse shock (with a rclativistic upstream to downstream 4- velocity, ur,s 3> 1, or 734 3> 1 
in the notation of ZK05) also outside its range of applicability, while the result for a fast-magnetosonic wave 
is approached in the opposite limit of a weak reverse shock (wrs C 1 or 734 - 1< 1). 
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shell and its deceleration due to the interaction with the e x terna l medium, following paper 
I and treating the latter part in more detail. iLevinsonl (120101 ) also identified regimes I 



and II. His expressions for the maximal Lorentz factor of the shell in regime II are only a 
factor of 1.09 lower than Eq. f lTTj) of the current work for k — 0, and a factor of 1.57 lower 
for k = 2 (the latter difference might appear larger since he used z = 0, E53 = 0.1 and 
A* m 33 for his fiducial values while the current work uses z = 2, E 53 = 1 and A* = 1). 
The current work finds that for k < 2 the maximal value of (r) is obtained at R u , and is a 
factor of ~ (<To/r cr )( 2-fc )/( 10 ~ 3fc ) > 1 (see Eq. 1491) larger than r cr . However, for the values 
of do < 10 3 and T cr = 180 considered by iLevinson (hoioh this factor if < 1.4 for k = 0, and 
thus consistent with his results (see his Fig. 5). 



Levinsonl ( 120101 ) has argued, however, that multiple sub-shells with an initial separation 



comparable to their initial width would collide and merge while still hig hly magnetize d, which 
is incorrect and in contradiction with paper I. An accompanying paper ( )Granotll201ll ) focuses 
on the possible role of multiple sub-shells, which can alleviate some of the requirements on 
the Lorentz factor of the outflow and may help accommodate GRB observations much better. 



6. Discussion and Conclusions 

This work has presented a detailed and unified treatment of the magnetic acceleration 
of an impulsive, initially highly magnetized (o"o ^> 1) shell and its deceleration by an unmag- 
netized external medium (with a power-law density profile). The dynamics divide into three 
main regimes (I, II, and III) and two more "exotic" regimes (relevant for an external density 
that either sharply drops with radius [II*], or is very large [IV], leading to a Newtonian flow). 

In regime I the external density is low enough that the shell accelerates almost as if into 
vacuum. At the coasting radius, R c ~ Rq(Tq, it reaches its maximal Lorentz factor of (T) ~ 
do < T C r (where T cr is given in Eq. [11]) and becomes kinetic energy dominate d. At R > R r 
this regime reverts back to the well-studied unmagnetized "thin shell" case (jSari &: Piran 



19951 ). where the shell coasts and spreads radially, A(R > R c ) ~ (R/R C )R , as its magne- 
tization rapidly decreases well below unity, a(R > R c ) ~ R c /R. In this regime the reverse 
shock is initially Newtonian, starts dominating the pressure behind the CD at i?Rs ~ R C r, 
and becomes mildly relativistic when it finishes crossing the shell, at i?r ~ Rdcc- The de- 
celeration radius, -Rdec, which corresponds to an observed deceleration time T dec , is where 
most of the energy dissipation in the shell takes place and most of the energy is transfered 
to the shocked external medium. Thus, both the reverse shock emission and the afterglow 
emission are expected to peak on the timescale of T ~ Td ec - At R > Rd cc (or T > Ta ec ) 
the flow quickly approaches the BM76 self-similar solution, which for GRBs signals the start 
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of the usual long-lived decaying afterglow emission. The magnetization at the deceleration 
radius is low, (a)(R dec ) ~ R c /R dcc ~ (cr /T cr ) 2 ( 4 ~ fc) /( 3 ~ fc) 1. If it is very low then mag- 
netic field amplification in the mildly relativistic collisionless (reverse) shock that develops 
could bring the downstream magnetic field to within a few percent of equipartition, thus 
allowing a good radiative efficiency for synchrotron emission, resulting in a bright reverse 
shock emission. In this regime the reverse shock emission and the afterglow emission both 
peak on a timescale T dec that is larger than the duration Tgrb of the prompt GRB emission, 
T dec /T GRB ~ Rdec/Rc 3> 1- Moreover, the degree of magnetization behind the reverse shock, 
~ (o")(-Rdcc) ~ ^GRB/Tdec <C 1, can be directly inferred from the ratio of these two observable 
times. 

In regimes II or III the shell remains highly magnetized near the deceleration radius, 
(cr)(i?dcc) ~ 00/Tcr ;§> 1, which strongly suppresses the reverse shock (which either becomes 
very weak or non-existent) and its associated emission. The energy in the flow is transfered to 
the shocked external medium (mostly near Rd ec ~ -Rcr) with very little dissipation within the 
original shell as long as ideal MHD holds. This is a highly magnetized "thick shell" case, and 
the afterglow onset time is similar to the initial shell light crossing time, T dec ~ (1 + z)Rq/c. 
For a single shell ejected from the source the prompt emission in this case might either arise 
from the onset of the forward shock emission (for an external shock origin, which makes it 
difficult to account for significant variability, and in which case Tgrb ~ Td ec ~ (1 + z)R /c) 
or alternatively due to magnetic reconnection events within the highly magnetized shell. The 
latter might be induced by the deceleration of the shell due to the external medium, in which 
case they might peak near R dec ~ R cr , since the angular size of causally connected regions 
(~ r _1 ) grows as the shell decelerates (Fcd oc _R( fc ~ 2 )/ 4 decreases with radius for k < 2) and 
at R < i?dec ~ -Rcr most of the energy is still in the original magnetized shell (this would 
again lead to T GRB ~ T dec ~ (1 + z)Rq/c). 

For the single shell case that was analyzed in this work there is either the low magnetiza- 
tion "thin shell" (regime I) or the high magnetization "thick shell" (regimes II or III) . There 
is no low magnetization "thick shell" case where a strong highly relativistic reverse shock 
develops, which can result in a bright reverse shock emission on a timescale comparable to 
that of the prompt gamma-ray emission in GRBs (T dec ~ Tgrb)- Similarly, there is no high 
magnetization "thin shell" case where the reverse shock is severely suppressed by a high 
magnetization in the shell near the deceleration radius and the afterglow onset occurs well 
after the prompt GRB emission (Td ec ^> Tgrb)- Therefore, a bright reverse shock emission is 
possible only in the low magnetization "thin shell" case - regime I, in which case this reverse 
shock emission (as well as the afterglow emission) would peak on a timescale larger than 
the du ration of the p rompt GRB emission, Td ec S> Tgrb ~ (1 + z)Rq/c. An accompanying 



paper ( jGranotll201ll ). however, shows that if the flow consists of many distinct sub-shells 
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instead of a single shell, then this may also allow a low magnetization "thick shell" regime. 



The Lorentz factor of the emitting region in GRBs must be high enough t o overcome the 
compactness problem and avoid excessive pair production within the sour c e (IKrolik fc Pier 



1991 



Fenimore et al. 1993; Woods fc Loeb 



1995 



Baring fc Harding! 119971 ; iLithwick fc Sari 



200ll ). It had been recently argued (lLevinsonll2010l ) that the interaction with the external 
medium might not enable an impulsive highly magnetized outflow in GRBs to accelerate up 
to sufficiently high Lorentz factors, and that its maximal achievable Lorentz factor is largely 
limited to T < r cr . This would pose a particularly severe problem for a stellar wind-like 
external medium (k = 2) for which typically r cr < 10 2 (see Eq. [H]). Recent high-energy 
observations by the Fermi Large Area Telescope (LAT) have set a lower limit of r > 10 3 for 



the em it ting region in a numbe r of GRBs with a bright high-energy emission (lAbdo et al. 



2009al lbl:l Ackermann et al.ll2010l ) using a simplified one- zone model. Howe yer, a more detailed 
and realistic treatment shows that the limit is lower by a factor of ~ 3 (IGranot et al.l 12008 



Ackermann et al.ll201ll ; lHascoet et al.ll201ll ). which would correspond to T > 10 2,5 for the 
brightest Fermi LAT GRBs. Nevertheless, this might still pose a problem for a single highly 
magnetized shell in a stellar- wind environment. 

Regime I implies a maximal Lorentz factor of the flow, r < oo <C r cr (where r cr is given 
in Eq. [H]). In regime II, a higher maximal Lorentz factor is possible for k < 2, and (r) 
peaks at T u where it exceeds r cr by a factor of ~ (cro/r cr )( 2_fc ^( 10_3fc ) > 1 (see Eq. [49]), while 
Tcd can reach values as high as ~ a at R < i? x (however, the material with such a Lorentz 
factor would carry only a small fraction of the total energy, ~ {Y CT / gq)^ 2 ^^ 2 ^^ < 1 at 
R ~ Ri). In regime III the typical Lorentz factor is close to that of the CD, (r) ~ Tcd, 

(2—k) II 

and for k < 2 they both peak at Rq where they exceed r cr by a factor of ~ Tc r > 1, 
while (r) ~ r cr {R/R cr Y k - 2 V* at Rq < R < R cr . All this could help increase (r) above T cr 
for k < 2. However, for a stellar wind-like (or steeper) external medium, k = 2 (or k < 2), 
we have (r) < T cr , which makes it very difficult to satisfy the observational constraints on T 
from pair opacity (mentioned above), and to a lesser ext ent also those from the onset time 
of th e afterglow emission (usually around a few hundred; ISari fc Piranlll999l : iNakar &; Piran 



2005 



Molinari et al.l 120071 ; IZou fc Piranl I201CH : iGruber et al.l |2011[ ) . The accompanying pa- 



per (lGranotll201ll ) shows, however, that if instead of a single shell the flow is initially divided 
into multiple, well separated sub-shells, then it can reach (r) ^> r cr and reasonably efficient 
internal shocks can naturally take place at such high Lorentz factors. This greatly helps to 
satisfy the observational constraints on T. 



The author thanks A. Spitkovsky, Y. E. Lyubarsky, T. Piran, A. Levinson and S. S. 
Komissarov for useful comments on the manuscript. This research was supported by the 
ERC advanced research grant "GRBs" . 
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Notation Definition Eq./Sect. 

To, Ao Initial Lorentz factor and lab-frame width of the unmagnetized shell §[2] 

pi = Ar~ k External medium rest mass density (r is the distance from the origin) §[2] 

p4 Proper rest mass density of the unmagnetized shell § [2] 

/ = pi/p\ Unmagnetized shell to external proper rest mass density ratio Eq. J2j 

Ig Sedov length (or radius) Eq. J3j 

Rm ~ min(i?r , Rn.o) ■ Radius where the reverse shock becomes Newtonian or relativistic §[2] 

i?p Radius where an external rest mass of E/TqC 2 is swept up Eq. JSJ 

i?A ~ max(ijp, R& q) . Radius where the reverse shock finishes crossing the unmagnetized shell §[2] 

R s ~ AqTq Radius where the shell starts spreading radially significantly §[2] 

RjVO> ^A,o Initial values (without radial spreading of the shell) of ijjy and R& Eq. J4]l, §[2] 

T = To(Ao/A) Reverse shock strength parameter (Newtonian for T > 1, rel. for T < 1) §[2] 

To Initial value (without radial spreading of the shell) of T Eqs. (5), (1 1 2 P 

r cr , A cr Critical Lorentz factor and width of the shell, respectively Eqs. HOP . lilt 

T, t Time when photons reach the observer and lab-frame time, respectively § El § El 

E, E cx t Total energy and energy in the shocked external medium, respectively § [2] § [3] § 

^BMi -f^EM 0: ^4dn ■ • ■ • Electromagnetic, initial electromagnetic and kinetic energies, respectively §E1 §H1 

<To = Bq/4-kpoc 2 ^> 1 . Initial value of the magnetization parameter §[3] 

B, Bo Lab-frame magnetic field and its initial value (at i?o), respectively §[3] §|4] 

po Initial proper rest mass density of the magnetized shell (at Ro) § [3j [4] 

/3ms, 0i Ums,o, r m s,o ■ • ■ Initial fast magnetosonic dimensionless speed, 4-velocity and Lorentz factor §[3] 

/3ms, »ms, r ms Fast magnetosonic dimensionless speed, 4-velocity and Lorentz factor §[3] 

r, a Lorentz factor and magnetization parameter of the shell, respectively § [3] 

(r), (a) Typical values of T and a - weighted means over the lab-frame energy §[3] 

i?0 ~ cto ~ Ao Initial radius (or lab-frame width) of the magnetized shell § [3] 

R c fa ct c ~ Rocr 2 , Coasting radius where the shell becomes kinetically dominated (in vacuum) §[3] 

Rco, Tcd Radius and Lorentz factor of the contact discontinuity (CD) that 

separates between the magnetized shell and the shocked external medium §[4] 

.Rshj r s h Radius and Lorentz factor of the shock front for the external shock §[4] 

§ = r/ct = x/ct Similarity variable §[4] Fig.[T] 

Ccd, £sh Values of § corresponding, respectively, to Rcd and R s h §[U Fig.[T] 

Value of J where the uniform region 3 in the Riemann problem starts § [4] Fig. [T] 

£ r f = — I3ms,0 Value of § at the head of the self-similar rarefaction wave §[3] Fig.[T] 

X, XCD Similarity variable of Blandford & McKcc (1976) and its value at Rev §H1 Eq. H20t 

Value of £ at the head of the secondary (or "reflected") rarefaction wave §[4] 

fit Dimensionless speed of the secondary ("reflected") rarefaction wave head Eq. 1221) 

a Ratio of pressure at the CD for the BM76 solution and a uniform region 2 Eqs. 123P - 12 i I 

R u tss ct u Radius where the secondary rarefaction wave reaches region 3, £»(tu) = £u Eqs. d 3 ll> . J40 D 

p = p/po = cr/o"0 Normalized shell proper rest mass density (or magnetization) §[4] 

/o = Po/ Pi{Ro) Initial (at Ro) magnetized shell to external proper rest mass density ratio §[4] 

Ri Radius where a = 1 just behind the contact discontinuity (CD) Eq. 1281 1 

o"u ~ crcDi pu = — ■ ■ Values of a and p, respectively, at £ = £ u (and also just behind the CD) Eqs. i l 2 r >|) . 1291 1. (130 P 

R cr ~ Ro^cr Critical radius where (r) reaches r cr in regimes II and III Eqs. 1461 1. 1471 1 

i? dcc Deceleration radius where most of the energy is transfered to the shocked 

external medium §[4] 

Rug Radius where a strong reverse shock develops in regime I Eq. 13411 

Rt.cu ~ ctcD Radius where the secondary rarefaction's head reaches the CD (£„ = £cd) Eqs. 1321 1. |43t 

Mrs Reverse shock upstream to downstream relative 4-velocity § 14.21 

L fs Ec/2Ro, ^cd • ■ • • Shell's mean total energy flux through a static sphere & its value at i?cD § 14.31 

r u = r C D{Ru) The CD as well as the typical Lorentz factor at Ru, F u ~ (T)(R U ) Eqs. (3DJ, H9t 

Tgrb = (1 + *)Ao/c. . Observed duration of the prompt GRB emission §[2] §|6] 

Tdcc ~ max (TGRB>Tr) Duration of peak reverse shock or afterglow emission (deceleration time) Eq. §|5] 



Table 1 : Notation and Definition of Some Quantities Used Throughout This Work. 



regime 


ordering of critical radii 


/o = po/pi(Ro) 




(TO 


I 

II 

III 

IV 


R < (Ri <) f R c < Ru ~ R*,cu < Rrs ~ Rcr < R doc ~ Rr 

R < (Rl <)* R u < Rcr ~ -Rdcc ~ R,,CD < Rc 

Ro ~ < Rcr ~ Rdcc ~ R*,CD < Rc 
Rdcc ~ Ro i *decAo ~ Tcr 4 ^> 1 


/o » a 7 " 2 * > 1 

r?r 2fc » f » r c 7 r - 2fc » i 

T < /() < cr () 

« A. « r c 7 r 2fc 

< /o < (To /3 

fo <r~, r cr > i 

fo « cr- 1 < 1 
/o « r?r 2fe « i 


r cr > o-„ > i 
/o 7 ^ » r cr » » i 

2 

^12-3 ^ ^ (Tq 

/cT^ « r cr « / A 
i « r cr « ct„ 12_3 ' c 

T cr »max(l,/ ( f r¥ ) 

r cr < i< (t 
/q — < rcr <<; 1 


r cr 

1 « ^0 « ff^ 

12- 3fc 

r cr < cto < r cr 2 

/O 7 ^ « (TO « fl 

12-3fc 

(to > r cr 2 > i 

(to > max(/,f,/ l 7 1 ) 

(T > i > r cr 

K (TO < / _1 



Table 2: The different regimes for k < 10/3 expressed in terms of fo = po/pi(Ro), r cr <~ (.foco) 1 ^ 8 2k ^ and ao- 
t This ordering of i?i is valid only for 2 < k < 10/3 and cr 7 " 2fe < / < erg <^=> r^~ fc)/2 < cr < T cr . 
* This ordering of #i is valid only for fc < 2 and erg < / < <r 7 Q ~ 2k T cr < er < lir~ fc)/2 . 



regime 


ordering of critical radii 


fo = po/pi(Ro) 


Tcr 




I 

II* 

f<k<l 
| <fc<4 

III 

f<k<l 
f <fc<4 

IV 


^ Ro ^ Rdcp < -Rc < -Rcr 
* ij < -Rl ~ i?dcp < Rc < Rcr 

Ro <C -Rit <C i?l ~ Rdcp < -Rc < Rcr 

-Ro < -Rcr ^ Rdcp ~ Rdec 
-Rdcc ~ Ro i tdccAo ~ Tcr 4 S> 1 


/o > o 1 / 3 » 1 
rfr 2fc » fo » r c 7~ » 1 

7-2fc „ f _~ 1/3 

rlr 2fc < /o « r c V 

T < /() < cr () 

/o «rJr 2fc , r cr > 1 

f < cr- 1 < 1 

fo <s rf r a < i 


2 

Tcr > cr () 12 - 3fc > 1 

/ A » r cr » f*** » 1 

2 

K O < Tcr < O-0 2 ~ 3fe 
1 « /(f^ « Tcr « /(P^ 

r cr >max(/ A ,/ ( ^) 

i«r CI « o-o 

r cr >max(l,/ 7 - 2fc ) 

fl~ 2k » r cr > 1 
r cr < i 

i 

/ ~ «fcr«l 


12-3fc 

i< o-o < r cr 2 
1 < o < / 3 

12-3fc 

K T cr 2 < OO < Tcr 

1 « / 3 « CO « /(P^ 
cr > max(/ s ,/ ~) 

o-o > r cr > i 

a„>max(/ ( 7 1 ,/ ( f^) 

i«/ ( r 1 «^««/«r 5¥ 

O-o > 1 > Tcr 

1 <S oo <§C /o 1 



Table 3: The different regimes for 10/3 < k < 4 expressed in terms of / = p /pi(R ), T cr <~ (fo&o) 1 ^ S 2k) and er - 
t This ordering holds for / > erg cr < iir^ 72 . 

* This ordering holds for cr 1/3 < / < <^=> lir~ fc)/2 < cr < r< r 2 ~ 3fc)/2 . 



-44- 





-i 
I 




u.y 




u.o 




0.7 


o 


0.6 


Q. 






0.5 


Q_ 






0.4 




0.3 




0.2 




0.1 








region 5 



-1.5 



-1 




region 2 



^ sh region 1 
1 1.5 



5 = X / Ct 



Fig. 1. — The self-similar structure when a cold magnetized shell initially at rest (occupying 
x < at t = 0) accelerates into an unmagnetized external medium (initially at rest and 
occupying x > at t = 0). For concreteness, I show the proper density normalized by its 
initial value in the magnetized shell (p ~ the density in region 5), for a = 6 and pi/po — 0.08. 
Such a self-similar solution in planar symmetry, where £ = x/ct, also corresponds to a 
solution in spherical symmetry, where £ = r/ct and (x,b,p) — > (r,rb,r 2 p) (see Eq. |19j). 
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Fig. 2. — Evolution of the typical Lorentz factor of the flow (where most of the energy 
resides), (r), as a function of radius R for k < 2 and for different values of the initial 
magnetization cr (and po oc l/o"o) and fixed values of the initial time or length scale (to ~ 
Rq/c or R ), energy (E ~ Lt « LR /c or L), and external density {k and A or pi(Ro) = 
ARy k ), which imply fixed r cr and R cr . In most cases of interest r cr ^> 1, so this is assumed 
to be the case here. The green and purple lines correspond to regimes I (1 < a < T cr ) and 
II (r cr < cr < Tlr 4 ~ k ^ 2 ), respectively. In regime III (<t > Tlr 4 ~ k ^ 2 ), (T)(R > Ro) becomes 
independent of cr and follows the thin solid red and blue lines. (The particular slopes in 
this figure are plotted for k = 0, but the general scalings are clearly indicated). 
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Fig. 3. — The same as Fig. [2] but with the addition of the Lorentz factor of the contact 
discontinuity, Tcd (dashed- dotted lines), until it becomes similar to the typical Lorentz factor, 
(r) (solid lines). The two remain similar up to the deceleration radius, after which Tcd starts 
falling behind (r) (at which stage only (r) is shown in the figure for clarity; dashed lines). 
(The particular slopes in this plot are drawn for k = 0, but the general scalings are clearly 
indicated). 
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Fig. 4. — The same as Fig. |3]but for 2 < k < 10/3. 
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Fig. 5. — Evolution of the typical magnetization (a) of the outflow as a function of radius 
R, corresponding to Fig. [2] (i.e. for k < 2, where each of the solid lines originating at 
R = Rq corresponds to a different value of cr ). The different regimes identified in the 
text are plotted using lines of different colors: regime I (1 < a < T CT ) in green, regime 
II (r cr < (To < r„ 4 k ^ 2 ) in purple, and regime III (cr > T^ 4 k ^ 2 ) in cyan. The lines 
corresponding to relevant critical radii (some of which depend on ctq) are also shown. 




Fig. 6.— The same as Fig. [5] but with the addition of the magnetization just behind 
the contact discontinuity, <jcd (dashed- dotted lines), until it becomes similar to the typical 
magnetization, (a). 




Fig. 7. — The same as Fig. [6] but for 2 < k < 10/3. 
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Fig. 8. — Evolution of the typical Lorentz factor (r) of the flow (where most of the energy 
resides; thick solid lines) and the Lorentz factor of the contact discontinuity, Tcd (dashed- 
dotted lines), as a function of radius R, for k < 2 and for different values of the external 
density normalization (/ = p /pi(R ) or pi(Ro) = ARq 1 * or A) and fixed values of all of 
the other model parameters (oo 3> 1, k, po, Ro, and therefore also E and L), which imply 
a constant R c ~ -Ro°"o an< ^ varying r cr ~ (/o^o) 1 ^ 8-2 ^ an d R C r ~ -Ror^ r . The purple, 
green and cyan lines correspond, respectively, to regimes I (r cr > cr > 1 or R cr > R c ), II 
(a 2/(12 - 3fc) < T cr < o- or a 2{3 - k)/{i - k) < RJR C < 1), and III (1 < T cr < <x 2/(12 - 3/c) or 
1 < R cv /Rq < (Tq kS> )- The borderlines between these regimes are indicated by thin dashed 
red lines. Within regime II, the thin dashed green line is the border between the regions 
with and without a break in r CD (i?) at Ri > R . (The particular slopes in this plot are 
drawn for k = 0, but the general scalings are clearly indicated). 
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Fig. 9. — The evolution of the magnetization a with radius for k < 2, similar to Fig. El but 
for different values of the external density normalization (/ = Po/pi(Ro) or Pi(Ro) — AR,Q k 
or A) and fixed values of all of the other model parameters (a ^> 1, k, po, Rq, and therefore 
also E and L), which imply a constant R c ~ Ro<Jq and varying r cr ~ (/oOo) 1 ^ 8-2 *^ and 
-R cr ~ -RcTcr. (similar to Fig. [8]). (The particular slopes in this plot are drawn for k = 1, but 
the general scalings are clearly indicated). 
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Fig. 10. — Phase space diagrams of the different dynamical regimes: in the fo-o~o plane 
for k < 10/3 (top left panel), r cr -<7 plane for k < 10/3 (top right panel), /o~r cr plane 
for k < 10/3 (bottom left panel), and in the fo-o~o plane for 10/3 < k < 4 (bottom right 
panel). Each regime is labeled and denoted by a different color, and the borders between the 
different regimes are indicated (by labeled thick black lines). 



